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This  report  summarizes  experimental  and  theoretical  studies  carried  out 
during  the  period  of  November  1985  to  November  1986  on  laser  cladding  using  a 
quarternary  mixture  of  nickel,  chromium,  aluminum  and  hafnium.  The  laser  clad 
region  has  been  studied  to  determine  the  distribution  of  undissolved  hafnium 
and  hafnium-rich  precipitates.  Also,  the  solid  solution  obtained  during  laser 
cladding  has  been  examined  to  determine  the  extension  of  solubility  of  hafnium 
in  nickel.  Because  of  inherent  rapid  cooling  during  laser  cladding  process, 
many  metastable  hafnium-rich  phases  are  formed.  These  phases  have  also  been 
studied  in  this  work. 

Mlcrostructural  evolution  of  Nl-Cr-Al-Hf  clad  on  AISI  1016  substrate  is 
accepted  for  publication  by  Acta  Metallurgica.  Similar  work  in  IN718  sub¬ 
strate  is  published  in  the  Proceedings  of  the  "High  Temperature  Coating" 
conference  held  at  TMS-AIME  1986  Fall  Meeting.  Another  paper  on  IN  718  with 
oxidation  data  is  submitted  to  the  Journal  of  High  Temperature  Technology. 

For  determining  the  composition  of  extended  solid  solution  formed  due  to 
rapid  cooling  in  laser  cladding  a  mathematical  model  is  presented.  This  model 
considers  diffusion  mechanism  for  mass  transport  in  a  one-dimensional  semi- 
infinite  molten  pool  of  the  cladding  material,  from  which  heat  is  removed  by 
conduction  through  a  one-dimensional  semi-infinite  solid  substrate.  The  rate 
of  solidification  has  been  obtained  by  modeling  the  cooling  process  as  a  com¬ 
posite  medium  heat  transfer  problem.  The  discontinuity  of  the  concentration 
field  has  been  simulated  using  a  nonequlllbrlum  partition  coefficient  and  then 
a  non-similar  exact  solution  for  the  mass  transport  equation  has  been  obtained 
using  a  set  of  similarity  variables  which  can  be  derived  using  Lie  group 
theory.  Two  papers  are  published  on  this  work;  one  will  appear  in  the  J.  of 
Metals,  February  1987  and  the  other  In  the  J.  Appl,  Physics,  April  1987  issue. 
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I.  INTRODUCTION 


A.  Literature  Review 

Structural  alloys  requiring  exposure  to  an  aggressive  atmosphere  at 
elevated  temperature  often  need  protective  coatings  for  oxidation  resistance 
to  achieve  enhanced  service  lifetime.  M+Cr+Al+RE  (M  =  Ni,  Fe,  Co;  RE  = 
Reactive  Elements)  systems  are  widely  used  for  such  coatings.  These  coatings 
tend  to  form  A1203  rich  scales.  The  A1 2O3  is  the  protective  coating  of  choice 
[1,2]  since  1)  diffusion  through  the  oxide  scale  is  very  slow,  2)  volatility 
is  limited,  3)  growth  kinetics  Is  slow,  and  4)  it  is  relatively  inert  in  a 
high  temperature  oxidation  environment. 

However,  the  adhesion  of  A1 2O3  to  the  substrate  is  poor.  Attempts  to 
Improve  the  adhesion  have  largely  centered  on  the  reactive  element  (RE) 
effect.  Y  or  Hf  are  mainly  used  to  Improve  the  adherence  of  A1203  scale  on 
the  M  Cr  A1  Y  or  MCrAlHf  alloys  and  coatings.  There  are  various  mechanisms 
proposed  by  various  authors  regarding  the  beneficial  effect  of  RE  addition  to 
MCrAl  alloys.  Some  of  these  are  listed  below. 

1.  The  formation  of  mechanical  "pegs"  which  anchor  the  scale  to  the  substrate 
metal  alloy  [1,3], 

2.  The  rare  earth  dlspersolds  provide  alternative  vacancy  coalescence  sites 
and  thus  prevent  the  vacancy  coalescence  at  the  scale-metal  Interface  [4], 

3.  Increased  bonding  forces  are  formed  between  the  scale  and  substrate  through 
preferential  segregation  of  the  reactive  elements  at  the  Interface  [3], 

4.  The  Interaction  with  trace  elements  such  as  sulfur  forming  refractory 
phases  and  thus  limiting  the  sulfur  segregation  at  the  scale-substrate 
Interface  [2], 

5.  The  reduction  of  scale  growth  processes  [5], 
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6.  The  enhancement  of  scale  plasticity  by  structure  modification  [61, 

7.  The  reduction  of  thermomechanical  differences  between  the  scale  and  the 
substrate  by  the  formation  of  graded  oxide  layers  [71,  and, 

8.  Refractory  metal  (e.g.  Re)  addition  to  nickel  superalloys  inhibits  the  y' 
coarsening  process  and  thus  results  in  higher  temperature  capabilities  [8]. 

Although  there  are  controversies  regarding  which  particular  mechanism  or 
groups  of  mechanisms  are  responsible  for  improved  high  temperature  properties, 
there  is  a  general  consensus  that  the  addition  of  reactive  metals  improves  the 
high  temperature  service  life  and  that  they  perform  best  when  they  are  in  a 
finely  dispersed  form. 

Using  conventional  techniques,  it  is  very  difficult  to  Incorporate  RE 
elements  in  finely  dispersed  form  above  1  wt.  percent,  which  is  the  solid 
solubility  limit  [91.  The  Inherent  rapid  heating  and  cooling  rate  In  laser 
processing  can  be  effectively  used  to  produce  an  alloy  with  extended  solid 
solution  and  uniform  distribution  of  phases  [ 10-13).  The  laser  cladding  tech¬ 
nique  was  used  to  produce  Nl-Fe-Cr-Al-Hf  alloy  with  high  Hf  contents  compared 
with  alloys  produced  by  conventional  techniques  such  as  plasma  spraying  of 
prealloyed  powders.  This  Is  a  technique  where  a  mechanical  mixture  of  powders 
of  choice  Is  delivered  Into  a  laser  melted  pool  and  solidified  rapidly 
[10,11,14].  The  Increase  In  solid  solubilities  of  the  alloying  elements  with 
RE  depends  upon  the  processing  conditions  for  laser  cladding.  The  main  vari¬ 
ables  during  laser  processing  are  laser  power,  size  and  shape  of  laser  beams, 
scan  velocity  and  the  chemistry  and  metallurgy  of  the  substrate.  The  overall 
composition  and  microstructure  of  the  laser  clad  materials  was  also  determined 
by  the  degree  of  mixing  and  cooling  rate.  For  determining  the  effects  of 
various  laser  processing  conditions  and  to  understand  the  process  physics  of 
laser  cladding  mathematical  modeling  was  also  done. 
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Ourlng  laser  cladding  the  alloy  Is  formed  at  the  surface  at  a  non- 
equilibrium  cooling  rate.  There  are  not  very  many  theoretical  studies 
published  in  the  literature  to  model  laser  surface  cladding  and  alloying  at  a 
nonequilibrium  cooling  rate.  A  three-dimensional  heat  transfer  model  for  CW 
laser  material  processing  was  developed  by  Mazumder  and  Steen  [ 15 j .  Con¬ 
vection  arising  due  to  surface  tension  in  a  laser  melted  pool  was  considered 
by  Chan,  et  al.  [16].  They  studied  Its  effect  on  surface  velocity,  surface 
temperature,  pool  shape,  and  the  cooling  rate.  The  distribution  of  matter  by 
diffusion  and  convection,  after  it  is  delivered  to  such  a  pool  was  examined  by 
Chande  and  Mazumder  [13].  There  has  also  been  some  efforts  to  predict  the 
formation  of  extended  solid  solution  under  -apld  cooling  using  thermodynamic 
variables  such  as  free  energy  and  chemical  potential.  The  thermodynamic 
approach  to  solidification  has  been  examined  very  well  by  Baker  and  Cahn 
[17].  Boettlnger  and  Perepezko  [18]  have  discussed  the  process  of  rapid 
solidification  from  the  point  of  view  of  thermodynamics.  Boettlnger  et  al. 
[19]  have  used  the  response  function  approach  of  Baker  and  Cahn  [17]  and 
stability  analysis  for  microsegregation-free  solidification.  More  about 
rapidly  solidified  materials  can  be  found  In  Refs.  [14,20].  The  mathematical 
model  of  the  present  study  Is  concerned  with  the  determination  of  composition 
of  the  extended  solid  solution  formed  due  to  rapid  cooling  in  laser  cladding 
using  space-time  variable  approach. 

B.  Objectives 

The  objectives  of  the  project  for  which  this  report  is  prepared  are  to 
develop  Nl-based  superalloys  of  Nl-Cr-Al-Hf  with  a  maximum  solid  solubility  of 
Hf  In  it  by  laser  cladding  technique  In  order  to  generate  stable  and  protec¬ 
tive  films  which  should  provide  Improved  oxidation  resistance.  It  is  also 


expected  that  these  alloys  with  extended  solid  solution  of  Hf  will  increase 
the  high  temperature  stability  by  retarding  the  y'  coarsening  kinetics  and 
provide  fine  precipitates  and  thus  improve  the  intermediate  and  high  temper¬ 
ature  strength. 

In  order  to  understand  the  high  temperature  service  behavior  of  these 
alloys,  it  is  necessary  to  know  about  the  crystal  structure,  stability  and 
chemical  composition  of  the  different  phases  developed  during  the  laser 
cladding  process  and  correlate  this  structure  with  the  laser  processing 
conditions.  Moreover,  to  the  best  of  the  authors'  knowledge  nobody  has 
reported  the  microstructure  evolution  of  the  rapidly  solidified  Ni-Fe-Cr-Al-Hf 
alloy  with  extended  solid  solution  of  Hf. 

Besides  these,  the  goal  of  the  present  work  is  to  develop  mathematical 
model  to  study  the  effects  of  various  laser  processing  conditions  on  the 
extension  of  solid  solubility  observed  during  laser  cladding.  These 
objectives  include  the  following: 

(I)  Effect  of 

(a)  Laser  power 

(b)  Laser  beam  diameter 

(c)  Laser-cladding  Interaction  time 

(d)  Cladding  thickness 

(e)  Cladding  powder  feed  rate 

(f)  Cooling  rate 

(g)  Speed  of  the  substrate  relative  to  the  laser  beam 

(h)  Convection  In  the  cladding  melt 

on  the  extension  of  solid  solubility  in  the  cladding  alloy. 

(II)  Development  of  nonequilibrium  phase  diagrams. 
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II.  EXPERIMENTAL  STUDIES 

A.  Experimental  Procedure 

Three  different  substrates  (AISI  1016,  Inconel  718  and  GTD111  a  modified 
version  of  Rene-80)  were  used  for  laser  cladding.  Ni,  Cr,  Al,  and  Hf  powders 
in  the  ratio  of  10:5:1:1  were  used  as  mixed  powder  feed  for  cladding.  The 
size  of  Ni,  Cr,  Al  powder  was  about  5  urn  in  diameter  whereas  Hf  powder  was  of 
44  pm  in  maximum  diameter.  The  average  chemical  analysis  of  the  laser  clad 

region  was  found  by  electron  probe  micro  analysis  to  be  70%  Ni,  10%  Cr,  10% 

Fe,  5%  Al ,  and  5%  Hf  (all  in  weight  percents).  The  cladding  treatments  were 
carried  out  using  an  AVC0  HPL  10  kW  CW  CO2  laser  with  F7  cassegrain  optics  (as 
shown  in  Fig.  1(a)).  The  laser  was  operated  in  the  TEM*q^  mode*.  The  beam 
focused  by  cassegrain  optics  was  reflected  downward  toward  the  substrate  by  a 
flat  mirror.  The  laser  was  operated  at  the  range  of  5  to  9  kW.  Specimens 

(approximately  6.4  mm  thick)  were  traversed  relative  to  the  laser  beam  at  a 

speed  of  6.35  to  25.4  mm  per  second.  The  powder  was  delivered  to  the  area  of 
interaction  by  a  pneumatic  powder  delivery  system  with  a  feed  screw.  The 
powder  flow  was  regulated  by  varying  the  speed  and  changing  the  size  of  the 
feed  screw.  The  flow  rate  for  the  present  study  was  approximately  0.2 
gm/second.  Argon  gas  with  a  flow  rate  of  7.7  gm/sec  was  used  to  maintain  a 
steady  powder  flow  through  the  copper  tubing  leading  to  the  substrate.  Powder 
was  placed  on  the  substrate  just  prior  to  its  contact  with  the  laser. 

A  shielding  of  argon  gas  with  a  shielding  box  was  used  to  minimize 
surface  contamination  during  laser  processing.  The  argon  gas  which  directed 
the  powder  flow  into  the  molten  pool  of  substrate  was  also  used  in  plasma 
suppression  at  the  laser-substrate  interaction  point.  The  function  of  the 


*0onut-shaped  laser  beam  with  a  Gaussian  power  distribution  in  the  outer 
ring  and  none  at  the  hole  of  the  donut. 
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Figure  1(a)  Schematic  diagram  of  laser  surface  cladding 

(b)  Three-dimensional  view  of  the  cladding  and  the  substrate 
(c-)  Geometric  configuration  used  In  the  present  model.  (The 
model  substrate  and  the  solidification  of  the  model  cladding 
have  been  shown  after  rotating  the  pool  ABCA  (see  Fig.  1(a))  by 
90  degrees  clockwise) 


shielding  box  was  to  provide  the  inert  environment  near  the  laser-substrate 
interaction  point  which  also  reduced  clad  porosity  [lOj. 

Microstructural  observation  of  clad  specimens  was  carried  out  by  Auger 
electron  spectroscopy,  optical,  and  scanning  electron  microscopy.  After 
mechanical  and  chemical  polishing  to  a  thickness  of  0.3  mm,  3  mm  discs  were 
punched  from  the  clad  material.  Specimens  for  the  TEM  observations  were 
prepared  by  the  jet  polishing  technique  with  an  electrolyte  of  200  c.c.  2- 
Butoxy-ethanol ,  400  c.c.  methanol,  50  c.c.  perchloric  acid  at  35  volts  with 
the  electrolyte  at  a  temperature  of  -15°C.  Samples  were  observed  using  a 
Phillips  EM  430  microscope  (equipped  with  an  EDAX)  operated  at  300  kV. 

B.  Results  and  Discussion 
Iron  Substrate 

Optical  microstructure  of  the  laser  clad  Ni-Cr-Al-Hf  alloy  has  been 
studied  in  this  work.  It  shows  uniform  distribution  and  formation  of  second 
phase  in  the  laser  clad  matrix  as  well  as  at  the  grain  boundaries.  At  low 
magnification  it  is  difficult  to  resolve  the  undissolved  Hf  particles.  The 
second  phase  appeared  to  be  spherical  in  nature.  The  concentration  of  each  of 
the  elements  (Ni,  Cr,  A1 ,  Hf,  Fe)  in  the  laser  clad  region  was  found  by  the 
electron  probe  microanalysis.  The  concentration  profile  of  each  alloying 
element  across  the  laser  clad/substrate  reg  on  was  also  determined  and  found 
to  be  uniform. 

In  order  to  see  the  presence  and  distribution  of  undissolved  Hf  in  the 
laser  clad  region.  Scanning  Electron  Microscopy  (SEM)  was  carried  out.  The 
SEM  micrograph  shows  undissolved  Hf  particles  dispersea  throughout  the  matrix 
and  also  few  at  grain  boundaries.  The  undissolved  Hf  particles  were  confirmed 
by  the  SEM  x-ray  microanalysis  and  this  is  discussed  in  the  next  section. 


To  see  the  detailed  internal  structure  of  the  laser  clad  region,  STEM 
investigation  was  carried  out.  A  general  survey  of  the  thin  foils  shows 
uniformly  distributed  undissolved  Hf  particles  of  1  to  9  pm  size,  and  Hf-rich 
precipitates  in  the  nickel  rich  matrix  containing  high  density  of  dislocations 
and  y 1  precipitates.  STEM  microanalysis  from  the  undissolved  Hf  revealed  only 
Hf  peaks. 

The  size  and  volume  of  the  y1  precipitates  in  y  matrix  depend  upon  the 
laser  processing  conditions  i.e.  cooling  rates.  There  is  linear  increase  in 
the  size  of  y1  precipitates  with  increasing  laser  power  for  a  constant 
traverse  speed.  Similarly,  the  volume  fraction  of  the  y'  precipitate 
increases  with  increase  in  laser  power.  The  volume  fraction  of  y'  precipitate 
is  higher  at  low  traverse  speed  as  compared  with  high  traverse  speed  for  a 
given  set  of  laser  power.  This  could  be  explained  on  the  basis  of  laser  beam 
interaction  time  with  the  alloy,  i.e.  rate  of  cooling  during  rapid 
solidification.  For  a  given  laser  beam  diameter,  there  would  be  less 
interaction  time  for  an  alloy  made  at  a  traverse  speed  of  14.8  mm/s  as 
compared  with  an  alloy  made  at  a  traverse  speed  of  6.25  mm/s.  Therefore,  the 
rate  of  cooling  would  be  faster  in  the  former  case  (14.8  mm/s)  as  compared 
with  the  latter  case  (6.35  mm/s). 

At  this  stage  It  is  difficult  to  comment  about  the  solid  solubility  of  Fe 
in  the  matrix  and  of  Hf  In  the  y1  phase  as  the  precipitates  were  very 
small.  But,  It  has  been  found  that  by  increasing  the  cooling  rates 
(particularly  In  laser  processing  technique),  there  was  an  Increase  in  solid 
solubility  of  the  alloying  additions  in  the  matrix  as  well  as  in  the  precipi¬ 
tates.  The  cooling  rates  attainable  by  this  process  are  known  to  enhance 
elemental  homogeneity  and  to  extend,  in  some  cases,  elemental  concentrations 
in  the  alloy  beyond  the  normal  saturation  limits  of  the  equilibrium  phase 


diagrams.  Furthermore,  under  normal  rapid  quenching  of  Ni-based  superalloys 
(e.g.,  Ni-Cr-Al),  the  y'  precipitation  could  only  be  seen  after  suitable  long 
aging  heat  treatments.  In  spite  of  a  high  cooling  rate  (-10  °C/s)  occurred 
during  laser  cladding  process,  y'  precipitates  have  been  observed  in  the  laser 
clad  Ni-Cr-Al-Hf  alloys.  It  clearly  indicates  that  the  addition  of  Hf  to  the 
Ni-Cr-Al  alloy  enhances  the  activation  energy  or  driving  force  for  the 
precipitation  of  y'  precipitates.  The  volume  fraction  of  y1  precipitation 
depends  upon  the  solid  solubility  of  Hf  in  the  matrix.  Thus,  the  initial  size 
and  volume  fraction  of  y'  precipitates  and  increase  in  solid  solubility  of  Fe 
and  Hf  have  an  influence  on  the  oxidation  resistance  properties.  On  the  basis 
of  STEM  microanalysis  of  the  matrix,  the  possible  composition  of  y'  precipi¬ 
tates  would  be  (Ni,  Cr,  Fe^Al.Hf). 

The  general  survey  of  the  thin  foil  from  the  laser  clad  samples  (under 
different  laser  processing  conditions)  shows  the  formation  of  many  unknown 
meta-stable  phases  with  different  microchemistry.  Since  the  equilibrium  phase 
diagram  for  this  quaternary  Ni-Cr-Al-Hf  is  not  available,  and  the  volume 
fraction  of  each  phase  developed  was  very  small,  there  is  difficulty  in 
identifying  each  of  the  phases. 

The  size  of  the  undissolved  Hf  particles  varied  from  0.5  pm  to  10  pm  in 
diameter  depending  on  the  traverse  speed  and  laser  power.  The  undissolved  Hf 
particle  size  increases  with  traverse  speed  and  laser  power.  It  is  expected 
that  the  undissolved  Hf  particle  size  will  be  higher  for  higher  speed  since 
Increased  speed  leads  to  decreased  Interaction  time.  Interaction  time  is 
critical  for  dissolution  of  high  melting  point  (2423°K)  material  such  as  Hf. 
The  Increase  of  undissolved  Hf  particle  diameter  with  the  increase  of  power  is 
not  obvious.  However,  this  can  be  explained  on  the  basis  of  specific  energy 
input  for  laser  cladding  [ 14 | .  For  constant  traverse  speed  and  powder  feed- 
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rate.  Increased  power  leads  to  slower  cooling  rate  which  results  in  reduced 
surface  tension  driven  convection  1 16 J .  Mass  transport  in  laser  surface 
modification  with  mixed  powder  feed  depends  primarily  on  convection  (13). 
Therefore  reduced  convection  will  cause  reduced  dissolution  leading  to  higher 
undissolved  Hf  particle  size.  Scatter  in  particle  size  distribution  data  is 
also  higher  for  higher  laser  power.  This  is  because  of  the  fact  that  higher 
power  Increases  the  width  of  the  molten  pool  [16,21]  resulting  in  a  wide  range 
of  cooling  rates  depending  on  the  position  in  the  pool.  This  variation  of 
cooling  rate  also  leads  to  variation  of  undissolved  Hf  particle  size  due  to 
variation  of  convection.  Similar  studies  have  been  carried  out  for  Inconel 
718  substrate.  For  details  and  micrographs  of  these  studies  see  the  Ref.  [22] 
for  Iron  substrate  and  Ref.  [31]  for  IN  718  substrate.  They  are  enclosed  on 
APPENOICES  I  and  II,  respectively. 
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III.  MATHEMATICAL  MODEL 

A.  Introduction 

The  physical  process  of  laser  cladding  involves  pouring  of  cladding 
material  powder  onto  the  substrate  to  be  clad  by  using  a  powder  delivery  sys¬ 
tem  and  simultaneously  melting  the  powder  using  a  laser  beam.  The  pool  of  the 
cladding  melt  which  is  formed  just  below  the  laser  beam  solidifies  by  losing 
heat  to  the  surrounding  air,  adjacent  cladding  and  the  solid  substrate  as  it 
moves  away  from  the  laser  beam  due  to  the  motion  of  the  substrate.  This  pro¬ 
duces  a  layer  of  cladding  on  the  solid  substrate  and  a  liquid  cladding  pool 
below  the  laser  beam  (see  Fig.  1(a)).  Figure  1(a)  shows  an  experimental  set-up 
for  laser  surface  cladding.  The  area  ABCA  shows  the  liquid  pool  of  cladding 
melt  which  is  in  contact  with  the  solid  substrate  along  AB.  It  is  also  in 
contact  along  BC  with  the  cladding  formed  on  the  substrate  and  the  portion  CA 
of  the  pool  Is  exposed  to  an  inert  atmosphere  to  avoid  oxidation  of  the 
cladding  material.  The  shape  of  the  strip  of  cladding  will  depend  on  the 
properties  and  the  temperature  of  the  cladding  melt  and  the  substrate, 
relative  motion  of  the  cladding  with  respect  to  the  substrate  and  the 
characteristics  of  the  cladding  powder  delivery  system.  In  the  present  model 
the  strip  of  cladding  BCDEB  (see  Fig.  1(a))  has  been  considered  to  have  semi- 
cylindrical  shape  (see  Fig.  1(b))  and  based  on  this  the  cladding  radius  (rc) 
is  determined  (see  Eq.  (18)).  This  paper  is  concerned  with  heat  and  mass 
transfer  In  the  cladding  melt  ABCA  to  model  extended  solid  solution  during 
laser  cladding.  This  liquid  pool  solidifies  by  conducting  away  heat  to  the 
substrate  and  to  the  solid  cladding  across  AB  and  BC,  respectively.  Also,  it 
loses  some  energy  to  the  ambient  Inert  gas  across  the  free  surface  CA.  But 
the  loss  of  energy  to  the  Inert  gas  across  the  surface  CA  will  be  smaller  than 
the  heat  loss  due  to  conduction  in  the  substrate  and  the  cladding  BCDEB. 
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However,  the  cladding  BCDEB  is  formed  from  the  liquid  cladding  melt  just 
solidified  and  hence  its  temperature  distribution  is  almost  uniform  and  much 
higher  than  the  temperature  of  the  substrate  which  is  immediately  below  the 
pool  ABCA.  So  the  heat  loss  across  BC  is  much  less  compared  to  that  across 
AB.  Therefore,  we  can  consider  that  the  pool  of  cladding  melt  (ABCA)  loses 
energy  only  in  one  direction  through  the  surface  at  AB  and  to  simplify  the 
heat  transfer  calculations  we  assume  that  the  pool  extends  up  to  infinity 
along  AB  as  well  as  BC.  Due  to  this  assumption  the  cladding  melt  solidifies 
along  AB  and  this  freezing  front  moves  upward  in  the  direction  of  BC.  The 
solid  substrate  has  also  been  taken  to  be  semi-infinite.  The  geometric 
configuration  of  this  pool  has  been  shown  in  Fig.  1(c)  by  rotating  the  pool  90 
degrees  in  the  clockwise  direction.  In  this  figure  the  freezing  front  has 
been  shown  to  be  planar.  This  is  true  when  a  pure  metal  solidifies.  For  a 
cladding  melt  of  an  alloy  system,  the  freezing  front  develops  curvature  due  to 
surface  tension  and  also  there  could  be  dendrite  and  cellular  growth  at  the 
solid-liquid  interface.  Moreover,  rapid  quenching  may  lead  to  growth  due  to 
spinodal  decomposition  and  nucleation  in  the  bulk  of  the  liquid  phase.  The 
stability  criterion  of  Mullins  and  Sekerka  [231  of  a  planar  interface  during 
solidification  of  a  dilute  binary  alloy  shows  that  the  planar  interface  will 
be  unstable  for  both  the  nickel-hafnium  and  the  nickel-aluminum  systems 
considered  in  this  study.  It  was  found  that  there  are  dendritic  growths  for 
both  cases,  although  the  dendrites  are  less  prominent  for  the  nickel-aluminum 
alloy  compared  to  those  of  nickel-hafnium  system.  Thus  the  solid-liquid 
Interface  for  nickel-aluminum  system  can  be  considered  planar.  Also  for  those 
systems  where  the  size  of  dendrites  is  extremely  small,  the  Interface  can  be 
taken  to  be  planar.  Nevertheless,  the  dendrites  will  affect  the  diffusion  of 
solute  atoms  at  the  freezing  front.  To  simplify  the  mass  transfer  analysis. 


the  freezing  front  has  been  assumed  to  be  planar  because  the  dendrite  size  is 
negligible  for  systems  such  as  nickel-aluminum.  Moreover,  the  trend  predicted 
by  the  model  is  more  important  than  the  exact  numerical  values  due  to  the 
paucity  of  high  temperature  materials  data  and  simplifying  assumptions. 
Besides  this  the  present  model  has  been  developed  based  on  a  few  more 
assumptions  listed  below. 

(1)  The  thermal  conductivity  and  the  thermal  diffusivity  for  a  mixture  is 
the  sum  of  the  volume-averaged  value  of  the  respective  transport 
properties  of  each  element  of  the  mixture, 

(2)  The  mass  diffusivity  of  each  element  in  the  liquid  phase  is  the  aver¬ 
age  value  of  self-diffusivlty  over  the  room  temperature  and  the  in¬ 

itial  temperature  with  modified  activation  energy  for  the  mixture, 

(3)  There  is  no  diffusion  of  mass  In  the  solid  phase, 

(4)  The  solute  segregated  at  the  solid-liquid  Interface  moves  in  the 

liquid  phase  by  diffusion  only, 

(5)  The  concentration  of  solute  in  the  liquid  alloy  is  equal  to  the 

nominal  composition  of  the  cladding  powder  mixture. 

(6)  Change  In  solute  concentration  in  the  liquid  phase  does  not  alter  its 
freezing  point  but  It  affects  the  heat  flux  across  the  solid-liquid 
interface,  and 

(7)  Only  50  percent  of  the  laser  energy  has  been  assumed  to  be  absorbed 
by  the  cladding  material.  Studies^  show  that  the  amount  of  laser 
energy  absorbed  by  different  materials  is  37-60  percent. 

Under  these  assumptions  It  Is  required  to  solve  the  one-dimensional  heat 
conduction  equation  to  obtain  the  speed  of  the  freezing  front  which  is  then 
utilized  to  solve  for  the  distribution  of  solute  in  the  liquid  phase.  In 
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laser  processing  of  material  usually  the  laser  beam  or  the  workpiece  moves  at 
a  certain  speed  and  when  the  laser  melted  pool  is  fully  developed  the  speed  of 
the  freezing  front  and  that  of  the  laser  beam  or  the  workpiece  are  identi¬ 
cal.  This  eliminates  the  need  for  solving  the  energy  transport  equation  to 
determine  the  speed  of  solidification  for  fully  developed  pool.  In  the  case 
of  laser  surface  alloying  the  solidification  progresses  with  segregation  of 
solute  atoms  from  the  solid  to  the  liquid  phase  and  this  alters  the  concentra¬ 
tion  of  solute  at  the  solid-liquid  interface.  This  will  affect  the  speed  of 
the  freezing  front  as  well  as  the  freezing  point  of  the  liquid.  To  take  the 
former  effect  into  account  the  energy  transport  equation  has  been  solved  while 
the  latter  one  has  been  Ignored  in  this  model. 
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B.  Problem  Formulations 

The  governing  equations  for  energy  transport  are: 
(i)  Region  I  (solid  substrate) 


a2T 


1  _  1 


aT 


1. 


ax 


=  at  ;  t 


0, 


<  x  <  0 


(ii)Region  2  (solidified  cladding) 


3  Tp  .  aT- 

~t  vir;t>0,0sx  - s(t) 

ax2  “2  3t 


(ill)Region  3  (liquid  cladding) 


3  T3  .  aT- 

~  Tt*  t  -  0*  S(t)  <  x  <  - 
ax2  a3  3t 


(1) 


(2) 


(3) 
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Here  T^,  T2.  and  T3  are  the  temperatures  and  a^,  c^,  and  a3  are  thermal 
diffuslvltles  of  regions  1,  2,  and  3,  respectively.  S(t)  is  the  location  of 

the  solid-liquid  interface  at  time  t.  If  k^,  k2,  and  k3  are  taken  to  be 

thermal  conductivities  of  the  regions  1,  2,  and  3,  respectively  and  if  the 
ambient  temperature,  the  freezing  point  of  the  cladding  melt  and  the  initial 
cladding  pool  mean  temperature  are  denoted  by  TQ,  Tm,  and  Tf,  respectively 

then  the  auxiliary  conditions  for  the  above  heat  transfer  problem  can  be 

written  as 


T1(— .  t) 


T 

0 


(4a) 


T 1 ( x  *0 )  -  Tq 

aT.  k9  aT, 

3X  ■  3x  at  x  ■  0 

T1(0,t)  =  T2(0,t) 

T2(x,t)  *  T3(x,t)  =  at  x  *  S(t) 

aTp  hs 

k2  ax~  '  k3  alT  *  oL  dt  at  x  *  S(t) 


(4b) 

(4c) 

(4d) 

(4e) 

(4f) 


where  p  Is  the  density  of  the  cladding  melt  and  L  its  latent  heat  of 
solidification.  Note  that  the  change  in  density  due  to  solidification  has 
been  neglected  here. 


T3(-.t)  »  Tf 


(4g) 


1 

! 

< 

t 

t 


T3(x,0)  =  Tf 


(4h) 


S(0)  *  0  (41) 

The  amount  of  solute  rejected  by  the  solidifying  cladding  and  Its  transport  In 
the  liquid  phase  are  governed  by  the  following  diffusion  equation. 


aC  y/o3  aC  n  a2C 

at  "  ax  ‘  D  ~ 2 

3t  /t  3X  ax^ 

where  c  Is  the  concentration  of  solute  in  the  liquid  phase,  D  is  the  mass 
diffusion  coefficient  and  Y  is  a  root  of  the  transcendental  Eq.  (12).  If  CQ 
is  the  Initial  concentration  of  solute  in  the  liquid  phase  and  kg  denotes  the 
equilibrium  partition  coefficient  at  the  solid-liquid  interface  then  the 
auxiliary  conditions  for  the  above  mass  transfer  problem  are 

(6a) 

C(-,t)  =  CQ  (6b) 

Y/c,1  a  r  y/03 

0(1  +  8)  g  =  ~~  (k.  -  l)c  at  x  -  0  (6c) 

/t  3X  /t  e 

in  a  coordinate  frame  of  reference  moving  with  the  solid-liquid  interface.  In 
the  expression  (6c),  b  *  x/D*,  where  x  is  the  interatomic  spacing  and  D*  is 
the  Inter-dlffuslvity.  This  expression  Is  obtained  by  incorporating  Into  the 
mass  balance  equation  at  the  solid-liquid  Interface  an  expression  for  non- 
equilibrium  partition  coefficient  derived  by  Aziz  [25). 


C(x,o)  =*  CQ 


I 
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C.  Method  of  Solution 


The  multi-region  heat  transfer  problem  can  be  solved  by  the  method  used 


by  Carslaw  and  Jaeger  [24].  The  temperature  distributions  in  various  regions 


are  given  by 


TiW> 


a,k,(T  -  T  ) 

T0  *  ^frTTrfTTTT  11  *  Erf<ni>1 


T2(n2) 


T  .  “l¥Tm-V 

o  *  Erf  (y) 


[1  +^Erf  (n-) 
al*2  i 


T3(n3)  - 


Tf  +  Erfc(y)  Erfc(n3* 


where 


*  x/2/a^t,  i  *  1,2,3 


and  the  location  and  the  speed  of  the  freezing  front  are  given  by 


S(t)  =  2y  (a^t) 


S(t)  =  Y/a3//t 


Here  k2  and  are  taken  to  be  equal  to  k3  and  a 3  respectively  and  y  is  ob¬ 


tained  from  the  following  transcendental  equation: 


Wo  .  Wf 

rTTFfTTT  T-riTTm 


_  a,  2 
0  -r-  Lyer 
*3 
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After  knowing  the  rate  of  formation  of  the  solid  phase  the  distribution  of 


solute  in  the  liquid  phase  can  be  obtained  by  solving  the  mass  diffusion  prob¬ 


lem  described  by  Eq.  (5)  and  the  associated  conditions  (6a-c).  An  exact  solu¬ 


tion  for  the  above  problem  can  be  obtained  by  using  similarity  variables  which 


are  determined  using  Lie  group  theory.  Oetails  of  this  theory  and  its  appli¬ 


cation  for  the  solution  of  differential  equations  can  be  found  [29,30],  Using 


this  group  theoretic  approach  the  similarity  variables  can  be  found  to  be  (See 


Refs.  [26,27,28]). 


— 1 —  and  ,  n  =  0,1,2,...,° 

✓t  n  An 


where  An's  are  constants  to  be  determined  later.  To  satisfy  the  boundary  con¬ 


dition  (6c)  using  these  new  variables  the  principle  of  superposition  will  be 


used  by  defining  the  concentration  C  as 


C  *  z  A  to  — 
n=0  n  n  (/t)n 


Using  (13)  and  (14)  the  Eq.  (5)  and  the  auxiliary  conditions  (6a-c)  can  be 


written  as 


20  -f  *  (y  +  2,/o3)  *  n  *n  -  0 


to  (y  ♦  ®)  *  C  6 
vn  '  o  no 


(16a) 


dton  y/<*3  (kg-l) 


*n  '  <*  -  sno> 


Y/a,  (k  -1)  n-1  . 

- L_ _L_  z  (-l)n  J  to1  (16b) 

u  J-0  J 


where  «n0  is  the  Kronecker  delta  and  A^s  of  expression  (14)  are  given  by 
An  =  (sr/a3)n,  n=0,  1,  2,  . .  ® 


For  this  problem  four  terms  of  the  series  (14)  were  computed  by  determin¬ 
ing  iiig,  4^,  and  <i>3  from  Eq.  (15).  These  solutions  are  (see  Ref. 
(26,27,28]). 


*0(r)  -  CQ  [1 


_ G  Erfc  (/2  r) _ 

(/8A)  Exp(-4Y2  a3/D)  +  G  Erfc(2YA^70T 


^(r) 


G*„<ro>  -(’-2  -  r02)/2 

rotG  ' 


4*2  (r)  •  - 


GU0(r0)  -  »!<■■„)) 

roG  +  <ro2  -  *> 


-(r2  -  r02)/2 
re 


*3(r)  *  - 


GU0(r0)  -  ♦1(r0)  *  »?(r0)| 
<r02  -  1)6  +  r0(3  -  r„2) 


(r2  -  1) 


-(r‘ 


r02)/2 


where 


G  =  Y/a3  (ke  -  1)  /2/D,  r  =  (y  +  2Y/a3)//?D  and  rQ  =  Y/2a3/0 


Using  these  expressions  for  4»g,  4>2  and  ii»3  the  concentration  of  solute  in 

the  liquid  phase  was  computed  using  the  Eq.  (14)  and  then  the  corresponding 
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concentration  in  the  solid  phase  (C$)  was  obtained  using  the  following 
expression  derived  by  Aziz  [25]  for  nonequilibrium  partition  coefficient. 


Ce  S(t) 6  +  k 

c1---!—  e  <m 

S(t)fl  +  1 

In  the  present  model  the  inter-diffusivity  D*  is  taken  to  be  equal  to  the 
mass  diffusivity  and  the  interatomic  spacing  x  is  taken  to  be  4A  which  is  the 
average  diameter  of  an  atom.  Other  physical  properties  for  the  alloy  were 
computed  at  its  nominal  composition  using  the  following  relationships: 


o  =  l/2:  W./o, 
j=l  J  3 


where  n  =  total  number  of  elements  present  in  the  liquid  pool  of  cladding 


L  *  2:  W.  L. 
j=l  J  3 

n  W.  _ 

k AT)  =  p  E  -J-  k  (T) 
1  j=l  °j  1J 


ki  =  ;  k. (T)dT/(Tf  -  To) 


cd(T)  -  r  W  C  (T) 
p  j=l  3  Pj 


cP-(  yndT/(Tm  -  v 


n  W,  _ 

T  — J-  n  . 


1.(1")  =  P  I  <*4  1 

1 


=  J  31(T)dT/(Tf  -  Tq) 


y.v.v-./v.v: 


and  the  diffusion  coefficient  for  the  j-th  solute  in  the  liquid  phase  is 


Tf  - 

°j  ■  {  °oj  e  RT  dT'<Tf  -  To> 


where  DQj  is  the  self-diffusion  coefficient  of  the  j-th  element  and  Q  is  the 
activation  energy  for  the  mixture  given  by 


Q  = 


In  all  these  expressions  for  thermophysical  properties  of  the  cladding  powder 
mixture,  the  index  i  =  1,2,3  refers  to  the  solid  substrate,  solidified 

cladding  and  the  liquid  cladding  region,  respectively,  and  j  =  1,2,3 . 

refers  to  the  j-th  element  of  the  cladding  powder  mixture.  Here,  Wj  and  wj 
are  the  weight  fraction  and  the  mole  fraction  of  the  j-th  element,  respec¬ 
tively.  k.j  and  a .  .  denote  respectively  the  thermal  conductivity  and  the 

thermal  diffusivity  of  the  j-th  element  in  the  i-th  region.  C  is  the 

Pj 

specific  heat  of  the  j-th  element.  The  thickness  of  the  cladding  can  be 
determined  by  taking  an  over-all  energy  balance  which  yields  the  following 
expression  for  the  radius  (rc)  of  the  semi-cyl indrical  strip  of  cladding: 


r  .  f  2fP _ 1 

C  '  ll(Tf  -  Tq)  C  +  L}  W 


(18) 


where  f  is  the  fraction  of  laser  energy  absorbed  by  the  cladding  material,  P 
is  the  power  of  the  laser  beam  incident  on  the  material,  and  v  is  the  speed  of 
the  workpiece  which  can  be  related  to  the  volumetric  rate  (q)  of  powder 
delivery  by 
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Q  =  \  (19) 

assuming  that  all  of  the  delivered  powder  form  cladding  on  the  substrate. 

0.  Results  and  Discussion 

The  results  of  the  above  model  was  compared  with  experimental  data. 
Laser  cladding  was  performed  on  nickel  substrate  with  two  different  cladding 
materials  of  nominal  composition  74%  Ni  and  26%  Hf  in  one  case  and  74%  Ni  and 
26%  A1  in  the  other  case.  The  solid  substrate  was  clad  with  these  cladding 
materials  which  take  the  shape  of  an  approximately  semi -cylindrical  strip  of 
metal  on  the  substrate  (see  Fig.  1(b)).  Electron  Probe  Micro  Analysis  of 
these  samples  shows  the  concentration  of  Hf  and  A1  in  the  matrix  of  Ni  in 
excess  of  what  is  predicted  by  the  equilibrium  phase  diagram.  The  present 
model  predicts  the  composition  of  the  extended  solid  solution  quite  well. 
These  results  have  been  presented  in  the  Table  I  for  Tf  =  2035°K  for  the  N i -Hf 
alloy  and  T^  =  1907°K  for  Ni-Al  alloy. 

Table  I  Comparison  of  Theoretical  Results  with  Experimental 
Results  of  Extended  Solid  Solution 

Laser  Laser  Beam  Speed  of  Nominal  Composition  Composition  in  the  solid 

Power  Diameter  the  Work-  of  the  cladding  solution  (wt.  %) 

(Kw)  (mm)  piece  material  Theoretical  Experimental 

(Inch/min)  (wt.  %)  Results  Results 

5  3  50  74  Ni,  26  Hf  3.05  Hf  3.58-6.50  Hf 

5  3  50  74  Ni,  26  A1  27.3  A1  22  A1 

Note  that  there  is  large  fluctuation  in  Hf  data  based  on  two  analytical 

techniques  used;  Electron  probe  microanalysis  and  x-ray  dispersive  analysis. 

Fine  distribution  of  Hf  particles  and  beam  diameter  of  two  Instruments  are 
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probably  the  source  of  fluctuation.  The  prediction  of  composition  of  the 
extended  solid  solution  based  on  the  present  model  is  within  30*  of  the 
experimental  results.  This  discrepancy  may  be  due  to  some  of  the  assumptions 
of  this  model.  For  example,  solidification  has  been  assumed  to  take  place  at 
a  constant  temperature  which  has  been  taken  to  be  the  freezing  point  of  the 
cladding  material  at  its  nominal  composition.  But  it  is  well  known  that  the 
change  in  composition  in  the  liquid  phase  will  alter  its  freezing  point.  Also 
the  present  model  utilizes  an  expression  (17)  for  nonequilibrium  partition 
coefficient  which  is  applicable  to  dilute  solution.  Apart  from  these  the 
presence  of  a  two-phase  zone  between  the  solidus  and  the  liquidus  lines  and 
the  surface  tension  driven  flow  causing  convection  in  the  liquid  pool  will 
affect  the  mixing  of  solute  in  the  liquid  phase  and  thus  alter  its  composition 
in  the  solid  phase.  Paucity  of  high  temperature  liquid  metal  data  also 
contributes  to  the  numerical  error.  However,  the  objective  of  this  study  is 
to  determine  the  trend  and  thus  understand  the  underlying  process  physics. 

The  above  model  was  used  to  study  the  effect  of  various  important  process 
parameters  such  as  the  cooling  rate,  cladding  powder  delivery  rate  on  the  com¬ 
position  of  hafnium  in  nickel  and  aluminum  matrices.  Details  of  these  studies 
can  be  found  in  (26,27,28],  Ref.  (26]  in  enclosed  on  APPENDIX  III. 
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IV.  CONCLUSIONS 


A.  Current  Works 

By  laser  cladding  process  uniform  and  fine  distribution  of  undissolved 
Hf,  Hf  rich  precipitates  and  matrix  with  extended  solid  solution  of  Hf  in  the 
laser  clad  region  could  be  achieved.  These  are  supposed  to  contribute  towards 
an  increase  in  the  oxidation  resistance  properties  of  N1  base  alloys  at  ele¬ 
vated  temperatures  due  to  various  mechanisms  as  discussed  earlier.  The  volume 
fraction  and  size  of  y'  precipitates  depend  upon  the  laser  processing  con¬ 
ditions.  Because  of  inherent  rapid  cooling  during  laser  cladding  process, 
many  meta-stable  Hf  rich  phases  are  formed.  The  formation  of  many  meta-stable 
Hf  rich  phases  would  also  contribute  toward  an  increase  in  the  oxidation 
resistance  properties.  The  overall  oxidation  resistance  properties  of  the 
laser  clad  alloy  produced  by  under-focused  laser  beam  was  found  to  be  better 
than  the  alloy  produced  by  over-focused  laser  beam  as  well  as  Inconel  718 
substrate.  The  present  results  suggest  that  the  reactive  elements  would 
change  the  mechanism  of  the  scale  formation  as  a  competitive  reaction  is 
taking  place  for  the  oxidation.  These  reactive  oxidized  particles  will  act  as 
sinks  for  excess  vacancies  thus  inhibiting  their  condensation  at  the 
scale/matrix  Interface  and  possibly  enhancing  scale  adhesion. 

The  mathematical  model  of  this  study  examines  the  extension  of  solid 
solubility  based  on  the  transport  of  energy  and  mass.  Solute  transport  has 
been  considered  to  take  place  only  in  the  liquid  phase  while  the  energy 
transport  has  been  considered  in  both  solid  and  liquid  phases.  Effect  of 
nonequilibrium  cooling  rate  on  solute  segregation  at  the  freezing  front  has 
been  taken  Into  account  by  considering  a  nonequilibrium  partition 
coefficient.  Using  this  the  mass  transfer  problem  has  been  solved 
analytically  in  a  frame  of  reference  moving  with  the  freezing  front  with  a 


velocity  determined  from  the  heat  transfer  problem.  These  mathematical 
solutions  have  been  utilized  to  study  the  effect  of  various  process  parameters 
on  the  concentration  of  solute  in  an  alloy.  The  cooling  rate  is  found  to  drop 
sharply  within  one  second  after  the  solidification  starts  and  it  is  also 
higher  for  lower  initial  pool  mean  temperature.  The  hafnium  concentration  in 
the  nickel-hafnium  alloy  is  found  to  be  higher  for  higher  cooling  rate.  Also, 
it  was  found  that  more  laser  energy  is  required  to  obtain  an  alloy  of  a  given 
composition  at  a  higher  powder  delivery  rate  than  that  for  a  lower  rate  of 
powder  delivery. 

8.  Plans  for  Future  Works 

A  parallel  laser  cladding  has  been  done  on  GTD  111  substrate  and  micro- 
structural  evolution  of  the  laser  clad  samples  are  under  progress.  Under  a 
laser  processing  conditions,  microstructural  changes  in  the  laser  clad  alloys 
made  on  three  different  substrates  would  be  compared  and  will  provide  an 
additional  information  about  the  effect  of  alloying  elements  on  the  micro¬ 
structure.  Similarly,  the  oxidation  properties  and  thermal  stability  of  the 
alloys  would  be  compared  and  optimum  condition  would  be  set  up  for  the  best 
oxidation  resistance  properties. 

We  will  also  carry  out  the  following  tasks  to  develop  mathematical  models 
which  will  allow  us  to  predict  the  composition  of  the  alloy  formed  due  to 
nonequilibrium  cooling  in  laser  cladding: 

(1)  To  develop  an  expression  for  nonequllibrium  partitioning  for  real 
solution. 

(2)  Develop  a  one-dimensional  model  for  nonequilibrium  phase  diagram  by  using 
temperature-dependent  equilibrium  partition  coefficient  In  the  expression 
of  nonequllibrium  partition  coefficient. 
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(3)  We  will  extend  the  above  models  to  two-  and  three-dimensional  cases  of 
heat  and  mass  transport  and  to  multi -component  alloy  systems. 
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ABSTRACT 

Alloys  and  coatings  for  alloys  for  improved  high  temperature  service  life 
under  aggressive  atmosphere  are  of  great  contemporary  interest.  There  is  a 
general  consensus  that  addition  of  reactive  elements  such  as  Hf  will  provide 
many  beneficial  effects  for  such  alloys.  The  laser  cladding  technique  was 

used  to  produce  Ni-Fe-Cr-Al-Hf  alloys  with  extended  solid  solution  of  Hf.  A 

10  kW  C02  laser  with  mixed  powder  feed  was  used  for  laser  cladding.  Optical, 

scanning  electron  (SEM)  and  scanning  transmission  electron  (STEM)  microscopy 
were  employed  for  microstructural  evolution  of  alloys  produced  during  laser 
cladding  processes.  The  electron  probe  microanalysis  and  the  auger  electron 
spectroscopy  were  also  used  for  micro-chemical  analysis  of  different  phases. 
Microstructural  studies  revealed  a  high  degree  of  grain  refinement, 
considerable  Increase  In  solubility  of  Hf  in  matrix  and  Hf  rich  precipitates 
and  new  metastable  phases.  This  paper  will  report  the  microstructural 
development  in  this  laser  clad  Ni-Fe-Cr-Al-Hf  alloy. 

INTRODUCTION 

Structural  alloys  requiring  exposure  to  an  aggressive  atmosphere  at 

elevated  temperature  often  need  protective  coatings  for  oxidation  resistance 
to  achieve  enhanced  service  lifetime.  MCrAl.R.E.  (M  =  Ni,  Fe,  Co)  systems  are 
widely  used  for  such  coatings.  These  coatings  tend  to  form  A^Oj  rich 
scales.  The  A1203  is  the  protective  coating  of  choice  (1,2)  since:  1)  dif¬ 
fusion  through  the  oxide  scale  Is  very  slow,  2)  volatility  Is  limited, 
3)  growth  kinetics  Is  slow,  and  4)  It  is  relatively  inert  In  a  high  temper¬ 
ature  oxidation  environment. 

However,  the  adhesion  of  AI2O3  to  the  substrate  is  poor.  Attempts  to 
improve  the  adhesion  have  largely  centered  on  the  reactive  element  (RE) 
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effect.  Y  or  Hf  are  mainly  used  to  improve  the  adherence  of  A1 2O3  scale  on 
the  M  Cr  A1  Y  or  MCrAIHf  alloys  and  coatings.  There  are  various  mechanisms 
proposed  by  various  authors  regarding  the  beneficial  effect  of  R.E.  addition 
to  MCrAl  alloys.  Some  of  these  are  listed  below. 

1.  The  formation  of  mechanical  "pegs"  which  anchor  the  scale  to  the  substrate 
metal  alloy  [1,3], 

2.  The  rare  earth  dispersoids  provide  alternative  vacancy  coalescence  sites 
and  thus  prevent  the  vacancy  coalescence  at  the  scale-metal  interface  [4], 

3.  Increased  bonding  forces  are  formed  between  the  scale  and  substrate  through 
preferential  segregation  of  the  reactive  elements  at  the  interface  [31, 

4.  The  interaction  with  trace  elements  such  as  sulfur  forming  refractory 
phases  and  thus  limiting  the  sulfur  segregation  at  the  scale-substrate 
interface  [2], 

5.  The  reduction  of  scale  growth  processes  [5], 

6.  The  enhancement  of  scale  plasticity  by  structure  modification  [6], 

7.  The  reduction  of  thermomechanical  differences  between  the  scale  and  the 
substrate  by  the  formation  of  graded  oxide  layers  [7],  and, 

8.  Refractory  metal  (e.g.  Re)  addition  to  nickel  superalloys  inhibits  the  y' 
coarsening  process  and  thus  results  in  higher  temperature  capabilities  [8[. 

Although  there  are  controversies  regarding  which  particular  mechanism  or 
groups  of  mechanisms  are  responsible  for  improved  high  temperature  properties, 
there  is  a  general  consensus  that  the  addition  of  reactive  metals  improves  the 
high  temperature  service  life  and  that  they  perform  best  when  they  are  in  a 
finely  dispersed  form. 


Using  conventional  techniques,  it  is  very  difficult  to  incorporate  RE 
elements  in  finely  dispersed  form  above  1  wt.  percent,  which  is  the  solid 
solubility  limit  (9).  The  inherent  rapid  heating  and  cooling  rate  in  laser 
processing  can  be  effectively  used  to  produce  an  alloy  with  extended  solid 
solution  and  uniform  distribution  of  phases  [ 10-13).  The  laser  cladding  tech¬ 
nique  was  used  to  produce  Ni-Fe-Cr-Al-Hf  alloy  with  high  Hf  contents  compared 
with  alloys  produced  by  conventional  techniques  such  as  plasma  spraying  of 
prealloyed  powders.  This  is  a  technique  where  a  mechanical  mixture  of  powders 
of  choice  is  delivered  into  a  laser  melted  pool  and  solidified  rapidly 
[10,11,141.  The  increase  in  solid  solubilities  of  the  alloying  elements  with 
RE  depends  upon  the  processing  conditions  for  laser  cladding.  The  main  vari¬ 
ables  during  laser  processing  are  laser  power,  size  and  shape  of  laser  beams, 
scan  velocity  and  the  chemistry  and  metallurgy  of  the  substrate.  The  overall 
composition  and  microstructure  of  the  laser  clad  materials  was  also  determined 
by  the  degree  of  mixing  and  cooling  rate.  For  the  optimum  laser  processing 
condition  and  possible  application,  a  scheme  is  needed  to  select  the  best 
combination  for  a  given  workpiece.  The  main  objective  is  to  develop  Ni-based 
superalloys  of  Ni-Cr-Al-Hf  with  a  maximum  solid  solubility  of  Hf  in  it  by 
laser  cladding  technique  in  order  to  generate  stable  and  protective  films 
which  should  provide  Improved  oxidation  resistance.  It  is  also  expected  that 
these  alloys  with  extended  solid  solution  of  Hf  will  increase  the  high  temper¬ 
ature  stability  by  retarding  the  •*'  coarsening  kinetics  and  provide  fine  pre¬ 
cipitates  and  thus  Improve  the  Intermediate  and  high  temperature  strength. 

In  order  to  understand  the  high  temperature  service  behavior  of  these 
alloys.  It  Is  necessary  to  know  about  the  crystal  structure,  stability  and 
chemical  composition  of  the  different  phases  developed  during  the  laser 
cladding  process  and  correlate  this  structure  with  the  laser  processing 


conditions.  Moreover,  to  the  best  of  the  knowledge  of  the  authors  nobody  has 
reported  the  microstructure  evolution  of  the  rapidly  solidified  Ni-Fe-Cr-Al-Hf 
alloy  with  extended  solid  solution  of  Hf. 


EXPERIMENTAL  PROCEDURE 

AISI  1016  steel  was  used  as  the  substrate  material.  Ni,  Cr,  Al,  and  Hf 
powders  in  the  ratio  of  10:5:1:1  were  used  as  mixed  powder  feed  for 
cladding.  Size  of  Ni,  Cr,  Al  powders  used  were  about  5  urn  in  diameter  whereas 
Hf  powder  was  of  44  pm  in  maximum  diameter.  The  average  chemical  analysis  of 
the  laser  clad  region  was  found  by  electron  probe  micro  analysis  to  be  70 %  Ni, 
10%  Cr,  10%  Fe,  5%  Al,  and  5%  Hf  (all  in  weight  percents).  The  cladding 


treatments  were  carried  out  using  an  AVCO  HPL  10  kW  CW  C0£  laser  with  F7 


Cassegrain  optics  (as  shown  in  Fig.  1).  The  laser  was  operated  in  the  TEM*0^ 
mode*.  The  beam  focused  by  Cassegrain  optics  was  reflected  downward  toward 
the  substrate  by  a  flat  mirror.  The  laser  was  operated  typically  at 
approximately  5,  6  and  7  kW.  Specimens  (approximately  6.4  mm  thick)  were 
traversed  relative  to  the  laser  beam  at  a  speed  of  approximately  6.35  and  14.8 
mm  per  second.  The  powder  was  delivered  to  the  area  of  interaction  by  a 
pneumatic  powder  delivery  system  with  a  feed  screw.  The  powder  flow  was 
regulated  by  varying  the  speed  and  changing  the  size  of  the  feed  screw.  The 
flow  rate  for  the  present  study  was  approximately  0.2  gm/second.  Argon  gas 
with  a  flow  rate  of  7.7  gm/sec  was  used  to  maintain  a  steady  powder  flow 
through  the  copper  tubing  leading  to  the  substrate.  Powder  was  placed  on  the 


substrate  just  prior  to  its  contact  with  the  laser. 


♦Donut-shaped  laser  beam  with  a  Gaussian  power  distribution  in  the  outer 
ring  and  none  at  the  hole  of  the  donut. 
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A  shielding  of  argon  gas  with  a  shielding  box  was  used  to  minimize 
surface  contamination  during  laser  processing.  The  argon  gas  which  directed 
the  powder  flow  into  the  molten  pool  of  substrate  was  also  used  in  plasma 
suppression  at  the  laser-substrate  interaction  point.  The  function  of  the 
shielding  box  was  to  provide  the  inert  environment  near  the  laser-substrate 
interaction  point  which  also  reduced  clad  porosity  [10]. 

Microstructural  observation  of  clad  specimens  was  carried  out  by  Auger 
electron  spectroscopy,  optical,  and  scanning  electron  microscopy.  After 
mechanical  and  chemical  polishing  to  a  thickness  of  0.3  mm,  3  mm  discs  were 
punched  from  the  clad  material.  Specimens  for  the  TEM  observations  were 
prepared  by  the  jet  polishing  technique  with  an  electrolyte  of  200  c.c.  2- 
Butoxy-ethanol ,  400  c.c.  methanol,  50  c.c.  perchloric  acid  at  35  volts  with 
the  electrolyte  at  a  temperature  of  -15°C.  Samples  were  observed  using  a 
Phillips  EM  430  microscope  (equipped  with  an  EDAX)  operated  at  300  kV. 


RESULTS  AND  DISCUSSION 


Optical  microstructure  of  the  laser  clad  Ni-Cr-Al-Hf  alloy  is  presented 
in  Fig.  2.  It  shows  uniform  distribution  and  formation  of  second  phase  in  the 
laser  clad  matrix  as  well  as  at  the  grain  boundaries.  At  low  magnification  it 
is  difficult  to  resolve  the  undissolved  Hf  particles.  The  second  phase 
appeared  to  be  spherical  in  nature.  The  concentration  of  each  of  the  elements 
(HI,  Cr,  Al,  Hf,  Fe)  in  the  laser  clad  region  was  found  by  the  electron  probe 
microanalysis.  Figure  3  represents  the  concentration  profile  of  each  alloying 
elements  across  the  laser  clad/substrate  region.  The  distribution  of  each 
alloying  elements  in  the  laser  clad  region  was  found  to  be  uniform. 

In  order  to  see  the  presence  and  distribution  of  undissolved  Hf  in  the 
laser  clad  region.  Scanning  Electron  Microscopy  (SEM)  was  carried  out.  The 
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SEM  micrograph  shows  undissolved  Hf  particles  dispersed  throughout  the  matrix 


and  also  few  at  grain  boundaries  as  indicated  by  arrows  in  Fig.  4.  The 


undissolved  Hf  particles  were  confirmed  by  the  SEM  x-ray  microanalysis  and 


this  is  discussed  in  the  next  section. 


To  see  the  detailed  internal  structure  of  the  laser  clad  region,  STEM 


investigation  was  carried  out.  The  general  survey  of  the  thin  foils  show  the 


fine  and  uniform  distribution  of  undissolved  Hf  particles  with  1  to  9  pm 


size,  and  Hf  rich  precipitates  in  the  nickel  rich  matrix  containing  high 


density  of  dislocations  and  y'  precipitates  (Fig.  5).  STEM  microanalysis  from 


the  undissolved  Hf  revealed  only  Hf  peaks,  Fig.  6.  The  detailed  micro- 


structural  analysis  would  be  discussed  below.  Figures  7(a)  and  7(b)  represent 


bright  field  (bf)  and  dark  field  (df)  electron  micrographs  of  the  laser  clad 


matrix  region.  The  dark  field  micrograph  from  one  of  the  (001)  reflections 


show  the  presence  of  ordered  precipitates  in  the  f.c.c.  matrix.  The 


precipitate  has  been  identified  as  y'  type  precipitate  which  is  an  ordered 


f.c.c.  phase.  The  corresponding  diffraction  pattern  from  the  matrix  is 


presented  in  Fig.  7(c).  Figure  7(c)  shows  the  superlattice  reflections  from 


the  precipitates  in  addition  to  the  major  reflections  from  the  gamma  matrix. 


The  local  zone  axis  of  both  gamma  (y)  and  gamma  prime  ( y ' )  grains  was  [100]. 


This  suggests  that  the  growth  of  the  y1  precipitate  was  along  the  cube 


direction  of  the  matrix,  i.e.  [100]  , // ( 100 ]  .  From  the  dark  field  image,  it 

y  y 


is  difficult  to  deduce  conclusively  the  shapes  as  the  precipitates  are  very 


small.  However,  it  is  believed  that  they  are  y'  precipitates  with  rectangular 


shape.  The  size  and  volume  of  the  y'  precipitates  depends  upon  the  laser 


processing  conditions  i.e.  cooling  rates,  which  will  be  discussed  later.  The 


composition  of  the  matrix  was  determined  by  EDX  and  is  displayed  in  Fig.  8. 


The  matrix  is  mainly  Ni  with  small  amounts  of  Cr,  Fe,  Al,  and  Hf  in  it.  From 


*  '  *  •/  %  v  *,  s  ■  *  % 


the  microchemical  analysis  it  is  possible  to  comment  on  the  substitution  of  Ni 
and  A1  by  Cr,  Fe,  and  Hf,  respectively.  The  Ni-Al-Hf  phase  diagram  is  the 
most  significant  as  it  supplies  the  only  information  available  on  the  substi¬ 
tution  behavior  of  Hf.  In  this  diagram,  the  y'  phase  field  is  extended 
towards  N^Hf,  indicating  that  Hf  substitutes  for  Al.  Rawlings  and  Bevan  [15] 
reported  that  Si,  V,  Mn,  Nb,  Hf,  and  Ta  alloying  additions  substitute  for  the 
aluminum;  Co  and  Cu  substitute  for  nickel  and  Cr,  Fe,  W,  and  Mo  substitute  for 
both  species  of  the  gamma  prime  (Ni^Al)  phase.  The  Mossbauer  investigation  by 
Nicholls  [16]  also  reported  an  interesting  result.  The  proportion  of  Fe  atoms 
on  aluminum  and  nickel  sites  was  found  to  be  a  function  of  Fe  concentration. 
At  low  concentrations  (2.5  at  %)  approximately  three-quarters  of  the  Fe  atoms 
were  on  aluminum  sites,  whereas  at  high  concentrations  (10  at  %Fe),  there  were 
approximately  equal  numbers  of  Fe  atoms  on  aluminum  and  nickel  sites.  The 
substitution  behavior  is  determined  by  electronic  rather  than  size  factor 
considerations  [15],  The  addition  of  these  alloying  element  has  also  an 
influence  on  the  mechanical  properties  of  y'  precipitate.  Considerable 
strengthening  is  obtained  in  the  y'  precipitate  when  the  alloying  addition 
substitute  for  aluminum  and  has  a  large  misfit  parameter  [17],  Thus  the 
addition  of  Hf  and  Fe  to  the  nickel-base  alloys,  i.e.,  Ni-Cr-Al,  would  modify 
both  mechanical  properties  as  well  as  oxidation  resistance  properties. 

In  order  to  determine  the  size  and  volume  fraction  of  the  y'  precipicates 
present  in  the  matrix,  series  of  dark  field  TEM  micrographs  were  taken  under 
strong  two-beam  condition  (operating  one  of  the  <100  reflections)  at  high 
magnifications.  In  order  to  minimize  the  error,  TEM  of  all  laser  processed 
samples  were  carried  out  at  the  same  high  magnifications  (130,000X  and 
160.000X).  The  volume  fraction  Vy  of  the  precipitates  have  been  determined  by 
quantitative  transmission  electron  microscopy  utilizing  approximate 
corrections  for  overlap  as  follows  [18,19]: 


where  the  foil  thickness  t  of  each  samples  were  estimated  and  the  apparent 
volume  fraction  and  the  average  precipitate  size  D  were  measured  in  each 
case. 

The  size  and  volume  of  the  y'  precipitate  in  y  matrix  depends  upon  the 
laser  processing  conditions.  Figure  9  shows  that  there  is  linear  increase  in 
the  size  of  y1  precipitates  with  increasing  laser  beam  power  for  a  constant 
traverse  speed.  The  size  of  the  y1  precipitate  is  comparatively  larger  at  low 
traverse  speed  (6.35  mm  per  sec)  as  compared  with  the  high  traverse  speed 
(14.8  mm  per  sec).  Similarly,  the  volume  fraction  of  the  y'  precipitate  also 
increases  with  the  function  of  laser  beam  power  (Fig.  10).  The  volume  frac¬ 
tion  of  y1  precipitate  is  higher  at  low  traverse  speed  as  compared  with  high 
traverse  speed  for  a  given  set  of  laser  beam  power.  This  could  be  explained 
on  the  basis  of  laser  beam  interaction  time  with  the  alloy,  i.e.  rate  of 
cooling  during  rapid  solidification.  At  a  given  laser  beam  power  (e.g.  6  kW), 
there  would  be  less  interaction  time  for  an  alloy  made  at  a  traverse  speed  of 
14.8  mm/s  as  compared  with  an  alloy  made  at  a  traverse  speed  of  6.25  mm/s. 
Therefore,  the  rate  of  cooling  would  be  faster  in  the  former  case  (e.g.  6  kW, 
14.8  mm/s)  as  compared  with  the  later  case  (e.g.  6  kW,  6.35  mm/s). 

At  this  stage  it  Is  difficult  to  comment  about  the  solid  solubility  of  Fe 
in  the  matrix  and  of  Hf  in  the  y'  phase  as  the  precipitates  were  very 
small.  But,  it  has  been  found  that  by  increasing  the  cooling  rates 
(particularly  In  laser  processing  technique),  there  was  an  increase  in  solid 
solubility  of  the  alloying  additions  In  the  matrix  as  well  as  in  the  precipi¬ 
tates.  The  cooling  rates  attainable  by  this  process  are  known  to  enhance 


elemental  homogeneity  and  to  extend,  in  some  cases,  elemental  concentrations 
beyond  normal  saturation  limits.  Furthermore,  under  the  normal  rapid 
quenching  of  Ni-based  superalloys  (e.g.,  Ni-Cr-Al),  the  y'  precipitation  could 
only  be  seen  after  suitable  long  aging  heat  treatments.  In  spite  of  such  a 

5 

high  cooling  rate  (~10  °C/s)  occurred  during  laser  melting  and  solidification 
process,  y'  precipitates  have  been  observed  in  the  laser  clad  Ni-Cr-Al-Hf 
alloys.  It  clearly  indicates  that  the  addition  of  Hf  to  the  Ni-Cr-Al  alloy 
enhances  the  activation  energy  or  driving  force  for  the  precipitation 
of  y'  precipitates.  The  volume  fraction  of  y‘  precipitation  depends  upon  the 
solid  solubility  of  Hf  in  the  matrix.  Thus,  the  initial  size  and  volume 
fraction  of  y1  precipitate  and  increase  in  solid  solubility  of  Fe  and  Hf  has 
an  influence  on  the  oxidation  resistance  properties.  On  the  basis  of  STEM 
microanalysis  of  the  matrix,  the  possible  composition  of  y‘  precipitate  would 
be  (Ni,  Cr,  Fe)3(Al,Hf). 

The  general  survey  of  the  thin  foil  from  the  laser  clad  samples  (under 
different  laser  processing  conditions)  shows  the  formation  of  many  unknown 
meta-stable  phases  with  different  microchemistry.  Since  the  equilibrium  phase 
diagram  for  this  quaternary  Ni-Cr-Al-Hf  is  not  available,  and  the  volume 
fraction  of  each  phase  developed  was  very  small,  there  is  difficulty  in 
identifying  each  of  the  phases. 

Figure  11  displays  an  electron  micrograph  of  the  laser  clad  region 
showing  the  blocky  type  of  precipitate  in  matrix.  Figure  12(b)  is  the  dark 
field  electron  micrograph  taken  from  one  of  the  precipitate  reflections.  The 
corresponding  selected  area  diffraction  pattern  and  indexing  (considering 
pseudo  cubic  crystal  structure)  is  presented  in  Figs.  12(c)  and  12(d),  respec¬ 
tively.  Because  of  inherent  rapid  solidification  during  laser  process,  meta¬ 
stable  phases  contained  distortions  and  strains,  therefore,  identification  of 
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phase  is  complex.  In  addition,  the  phase  diagram  of  this  complex  alloy  system 
is  not  available  and  volume  fraction  of  each  phase  is  very  small. 

The  corresponding  x-ray  micro  analysis  from  the  precipitate  is  presented 
in  Fig.  12.  It  revealed  that  the  precipitate  is  rich  in  Hf  with  relatively 
small  amounts  of  Ni,  Cr  and  Fe  in  it.  As  the  precipitate  contains  mainly  Hf 
and  Ni  with  very  small  amount  of  Cr,  Fe  and  A1 ,  thus  pseudo  binary  N i -Hf 
diagram  can  be  considered.  Based  on  the  diffraction  analysis  and  micro¬ 
chemistry  of  the  precipitate,  it  appeared  that  the  precipitate  is  having 
tetragonal  crystal  structure.  Detailed  convergent-beam  microdiffraction 
diffraction  analysis  of  this  precipitate  will  be  reported  in  the  future. 

The  size  of  the  undissolved  Hf  particles  varied  from  0.5  ym  to  10  ym  in 
diameter  depending  on  the  traverse  speed  and  laser  power  as  shown  in  Fig.  13. 
Figure  13  shows  that  undissolved  Hf  particle  size  increases  with  traverse 
speed  and  laser  power.  It  is  expected  that  the  undissolved  Hf  particle  size 
will  be  higher  for  higher  speed  since  increased  speed  leads  to  decreased 
interaction  time.  Interaction  time  is  critical  for  dissolution  of  high 
melting  point  (2423°K)  material  such  as  Hf.  The  increase  of  undissolved  Hf 
particle  diameter  with  the  increase  of  power  is  not  obvious.  However,  this 
can  be  explained  on  the  basis  of  specific  energy  input  for  laser  cladding 
[14].  For  constant  traverse  speed  and  powder  feedrate,  increased  power  leads 
to  slower  cooling  rate  which  results  in  reduced  surface  tension  driven 
convection  [20].  Mass  transport  in  laser  surface  modification  with  mixed 
powder  feed  depends  primarily  on  convection  [13].  Therefore  reduced 
convection  will  cause  reduced  dissolution  leading  to  higher  undissolved  Hf 
particle  size.  Scatter  in  particle  size  distribution  data  is  also  higher  for 
higher  laser  power.  This  is  because  of  the  fact  that  higher  power  increases 
the  width  of  the  molten  pool  [20,21]  resulting  in  a  wide  range  of  cooling 
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rates  depending  on  the  position  in  the  pool.  This  variation  of  cooling  rate 
also  leads  to  variation  of  undissolved  Hf  particle  size  due  to  variation  of 
convection. 

It  would  be  appropriate  to  comment  here  that  by  increasing  the  RE 
addition  to  the  Ni  based  supper  alloys  it  is  expected  that  it  would  contribute 
to  the  improvement  of  the  oxidation  resistance  properties  following  many  of 
the  mechanisms  as  described  earlier. 
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SUMMARY 


By  laser  cladding  process  uniform  and  fine  distribution  of  undissolved 
Hf,  Hf  rich  precipitates  and  matrix  with  extended  solid  solution  of  Hf  in  the 
laser  clad  region  could  be  achieved.  These  are  supposed  to  contribute  towards 
an  increase  in  the  oxidation  resistance  properties  of  Ni  base  alloys  at 
elevated  temperatures  due  to  various  mechanisms  as  discussed  earlier.  The 
volume  fraction  and  size  of  y 1  precipitates  depend  upon  the  laser  processing 
conditions.  Because  of  inherent  rapid  cooling  during  laser  cladding  process, 
many  meta-stable  Hf  rich  phases  are  formed.  The  formation  of  many  meta-stable 
Hf  rich  phases  would  also  contribute  toward  an  increase  in  the  oxidation 
resistance  properties.  Thus,  overall  oxidation  resistance  properties  of  the 
alloy  developed  by  the  laser  cladding  process  would  be  expected  to  be  better 
than  Ni  base  alloy  made  by  the  conventional  technique. 
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Fig.  1  Optics  used  for  laser  cladding. 
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Fig.  3 


Electron  micro  probe  analysis  of  laser  clad  region  and  substrate 
showing  the  uniform  distribution  of  alloying  elements. 
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Fig.  6 


STEM  X-ray  microanalysis  from  undissolved  Hf 
particles  showing  mainly  Hf  peaks. 


.  8  STEM  X-ray  microanalysis  from  the  matrix  showing  that 
matrix  is  rich  In  HI  with  small  amount  of  Cr,  Fe,  A1 
and  Hf  alloying  elements  In  it. 
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Traverse  Speed  6.35  mm/sec 
Traverse  Speed  14.8  mm/sec 
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Traverse  Speed  6.35  mm/sec  O 
Traverse  Speed  14.8  mm/sec  A 


-►*  Laser  Power  (kW) 


10  A  plot  of  volume  fraction  of  y'  precipitate  with 
the  function  of  laser  power  and  traverse  speed. 
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12  STEM  X-ray  microanalysis  of  the  precipitate 

indicating  that  the  precipitate  is  Hf  rich  phase 
with  small  amount  of  Hi,  Cr,  and  Fe  in  it. 
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Abstract 


Insitu  Ni-Cr-AI-Hf  alloy  has  been  developed  by  laser  surface  cladding 
with  mixed  powder  feed  for  improved  high  temperature  oxidation  resistance. 
The  oxidation  resistant  materials  for  service  at  elevated  temperatures  must 
satisfy  two  requirements:  diffusion  through  oxide  scale  must  occur  at  the 
lowest  possible  rate  and  oxide  scale  must  resist  spallation.  The  formation  of 
an  A1 2O3  protective  scale  fulfills  the  former  requirements  but  its  adherence 
is  poor.  Reactive  metal  such  as  Hf  is  added  to  improve  adhesion.  A  10  kw  CO2 
laser  was  used  for  laser  cladding.  Optical,  SEM  and  STEM  microanalysis  tech¬ 
niques  were  employed  to  characterize  the  different  phases  produced  during 
laser  cladding  process.  Microstructural  studies  showed  a  high  degree  of  grain 
refinement,  increased  solid  solubility  of  Hf  in  matrix  and  Hf  rich  precipi¬ 
tates.  Thermo  gravimetric  analysis  was  carried  out  to  determine  the  oxidation 
properties  of  these  clad  alloys  with  extended  solid  solution  of  Hf.  Consider¬ 
able  improvement  over  the  base  metal  was  observed.  This  paper  discusses  the 
microstructural  development  in  this  laser  clad  alloy  and  its  effect  on 
oxidation. 

Keywords:  laser  cladding,  oxidation  properties,  high  temperature  coating. 


Introduction 


Alloy  selection  for  gas  turbine  materials  application  is  to  ensure 
superior  mechanical  properties  and  good  corrosion  and  oxidation  properties  at 
elevated  temperatures.  In  recent  years,  attention  has  been  focused  on 
developing  suitable  coatings  to  enhance  service  lives  and  are  applied  over 
less-oxidation  resistant  but  mechanically  stronger  alloy  substrates.  MCrAl 
(M  =  Ni,  Fe,  Co)  systems  are  widely  used  for  such  coatings.  These  coatings 
tend  to  form  AI2O3  rich  scale  which  acts  as  a  protective  coating  because: 

1)  diffusion  of  oxygen  through  the  oxide  scale  is  very  slow, 

2)  volatility  is  limited, 

3)  growth  kinetics  is  slow,  and 

4)  it  is  relatively  inert  in  a  high  temperature  oxidation  environment. 

It  has  been  well  established  that  the  addition  of  trace  level  of  reactive 
elements  such  as  yttrium,  hafnium,  etc.  to  coating  alloy  composition  greatly 
enhances  oxide  scale  adherence  [1-3], 

There  are  various  mechanisms  proposed  by  various  authors  regarding  the 
beneficial  effect  of  R.E.  addition  to  MCrAl  alloys.  Some  of  these  are  listed 
below. 

1.  The  formation  of  mechanical  "pegs"  which  anchor  the  scale  to  the  substrate 


metal  alloy  [1,31, 


2.  The  rare  earth  dispersoids  provide  alternative  vacancy  coalescence  sites 
and  thus  prevent  the  vacancy  coalescence  at  the  scale-metal  interface  [4], 

3.  Increased  bonding  forces  are  formed  between  the  scale  and  substrate  through 
preferential  segregation  of  the  reactive  elements  at  the  interface  [3], 

4.  The  interaction  with  trace  elements  such  as  sulfur  forming  refractory 
phases  and  thus  limiting  the  sulfur  segregation  at  the  scale-substrate 
interface  [5], 

5.  The  reduction  of  scale  growth  processes  [61, 

6.  The  enhancement  of  scale  plasticity  by  structure  modification  [7J, 

7.  The  reduction  of  thermomechanical  differences  between  the  scale  and  the 
substrate  by  the  formation  of  graded  oxide  layers  [8],  and, 

8.  Refractory  metal  (e.g.  Re)  addition  to  nickel  superalloys  inhibits  the  y' 
coarsening  process  and  thus  results  in  higher  temperature  capabilities  [9]. 

Using  conventional  techniques,  it  is  very  difficult  to  incorporate  rare 
earth  (RE)  elements  in  finely  dispersed  form  above  1  wtX,  which  is  the  solid 
solubility  limit  [10].  The  inherent  rapid  heating  and  cooling  rates  in  laser 
processing  can  be  effectively  used  to  produce  an  alloy  with  extended  solid 
solution  and  uniform  distribution  of  phases  [11-14].  The  laser  cladding  tech¬ 
nique  was  employed  to  produce  Nl-Cr-Al-Hf  alloy  with  high  Hf  contents  compared 
with  alloys  produced  by  conventional  techniques.  The  over  all  composition  and 
microstructure  of  the  laser  clad  materials  was  determined  by  the  degree  of 
mixing,  cooling  rate,  and  other  laser  processing  variables  [11]. 

The  objective  of  this  study  was  to  develop  an  alloys  of  Ni-Cr-Al-Hf  with 
increased  Hf  content  which  can  provide  better  high  temperature  oxidation 
resistance  properties.  To  the  authors  best  knowledge,  no  one  has  reported  the 
development  of  such  coatings  by  laser  cladding  process. 
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Experimental  Procedure 

The  main  objective  of  the  laser  cladding  process  is  to  clad  a  substrate 
with  another  material  having  a  different  chemistry  by  melting  a  thin  inter¬ 
facial  layer  to  produce  a  metallurgical  bond  with  minimum  dilution  of  the  clad 
layer  by  the  substrate  to  provide  improved  surface  properties.  The  clad 
material  can  be  added  either  by  predeposition  techniques  (e.g.  preplacing  a 
powder  on  the  substrate,  plasma  spraycoat,  electroplating)  or  co-deposition 
techniques  (e.g.  pneumatic  delivery  of  powder  into  the  laser  melted  pool,  wire 
feed  and  ribbon  feed).  In  this  study  a  mixture  of  powders  were  pneumatically 
delivered  into  the  laser  melted  pool.  This  technique  is  more  energy  efficient 
and  provides  higher  cooling  rate  (15]  which  is  needed  to  achieve  extended 
solid  solution. 

Inconel  718  plate  was  used  as  a  substrate  material.  Ni,  Cr,  Al,  and  Hf 
powders  in  the  ratio  of  10:5:1:1  were  used  as  mixed  powder  feed  for  clad¬ 
ding.  The  cladding  treatments  were  carried  out  using  an  AVCO  HPL  10  kW  CW  C02 
laser  with  F7  Cassegrain  optics  (as  shown  in  Fig.  1).  The  laser  was  operated 
in  the  TEM*q^  mode*.  The  beam  focused  by  Cassegrain  optics  was  reflected 
downward  toward  the  substrate  by  a  flat  mirror.  The  laser  was  operated 
typically  at  approximately  5,  6,  and  7  kW.  Specimens  (approximately  6.0  mm 
thick)  were  traversed  relative  to  the  laser  beam  at  a  speed  of  approximately 
6.35  and  10.6  mm  per  second.  The  powder  was  delivered  to  the  area  of  inter¬ 
action  by  a  pneumatic  powder  delivery  system  with  a  feed  screw.  The  powder 
flow  was  regulated  by  varying  the  speed  and  changing  the  size  of  the  feed 
screw.  The  flow  rate  for  the  present  study  was  approximately  0.2  gm/second. 


v  *Donut-shaped  laser  beam  with  a  Gaussian  power  distribution  in  the  outer  ring 

and  none  at  the  hole  of  the  donut. 
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Argon  gas  with  a  flow  rate  of  0.017  lb./s  was  used  to  maintain  a  steady  powder 
flow  through  the  copper  tubing  leading  to  the  substrate.  Powder  was  placed  on 
the  substrate  just  prior  to  its  contact  with  the  laser. 

A  shielding  of  argon  gas  with  a  shielding  box  was  used  to  minimize 
surface  contamination  during  laser  processing.  The  argon  gas  which  directed 
the  powder  flow  into  the  molten  pool  of  substrate  was  also  used  in  plasma 
suppression  at  the  laser-substrate  interaction  point.  The  function  of  the 
shielding  box  was  to  provide  the  inert  environment  near  the  laser-substrate 
interaction  point  which  also  reduced  clad  porosity. 

Microstructural  observation  of  clad  specimens  was  carried  out  by  optical 
and  transmission  electron  microscopy.  After  mechanical  and  chemical  polishing 
to  a  thickness  of  0.3  mm,  3  mm  discs  were  punched  from  the  clad  material. 
Specimens  for  the  TEM  observations  were  prepared  by  the  jet  polishing  tech¬ 
nique  with  an  electrolyte  of  200  c.c.  2-Butoxy-ethanol ,  400  c.c.  methanol,  50 
c.c.  perchloric  acid  at  35  volts  with  -15°C  temperature  of  electrolyte. 
Samples  were  observed  using  a  Phillips  430  microscope  (attached  with  an  EDAX) 
operated  at  300  kV. 

Oxidation  was  performed  on  a  Thermo  Gravimetric  Analyzer  (TGA)  of  Dupont 
1090  system.  Test  was  carried  out  up  to  temperature  1050°C  with  air  flow  rate 
between  45-50  cc/min  for  a  duration  of  60  minutes. 

Results 

Figure  2  represents  a  plot  of  weight  gain  with  function  of  time  for  a 
period  of  60  minutes  at  three  temperatures  (700,  900,  and  1050°C)  of  the  laser 
clad  samples  (laser  power  5  kW  and  traverse  speed  6.35  mm/sec).  It  revealed 
that  there  is  parabolic  weight  gain  of  the  sample  at  oxidation  temperature 
700°C,  and  then  it  decreases  at  oxidation  temperature  900°C.  There  was  not 


any  significant  change  in  the  weight  of  sample  observed  at  high  oxidation 
temperature  1050°C.  Figure  3  represents  the  oxidation  data  of  the  laser  clad 
sample  operated  at  different  laser  power  and  traverse  speed.  It  shows  that 
the  maximum  weight  gain  was  observed  in  the  laser  clad  sample  with  5  kW  laser 
power  and  6.35  mm/sec  traverse  speed  (with  over  focused  laser  beam)  compared 
with  the  other  samples  of  5,  6,  and  7  kW  laser  power,  and  traverse  speed  of 
6.35  and  10.58  mm/sec  (with  under  focused  laser  beam).  Figure  3  also  revealed 
that  keeping  the  traverse  speed  constant  (i.e.,  6.35  mm/sec),  with  increasing 
laser  power  for  under  focused  laser  beam  results  decrease  in  weight  gain. 
Table  I  gives  the  oxidation  data  for  Inconel  718  at  oxidation  temperature 
900°C  for  1  hour.  On  comparing  the  oxidation  data,  the  laser  clad  alloys  show 
better  oxidation  resistance  properties.  In  addition,  the  alloys  prepared  by 
under  focused  laser  beam  shows  a  dramatic  increase  in  the  oxidation  resistance 
properties  as  compared  with  the  substrate  and  alloys  prepared  by  over  focused 
laser  beam.  This  will  be  discussed  more  in  the  next  section. 

The  average  composition  of  each  elements  was  found  to  be  uniform  through¬ 
out  the  laser  clad  region  except  Hf.  Hf  distribution  was  found  to  be  uniform 
through  the  laser  clad  matrix  with  compositional  fluctuations.  This  depends 
upon  the  location  of  the  microprobe  analysis  on  either  of  Hf  rich  precipi¬ 
tates,  Hf  in  the  matrix  or  undissolved  Hf  particles  in  the  matrix,  thus  the 
composition  of  the  other  alloying  elements  changes  accordingly.  This  would  be 
discussed  in  more  detail  in  the  next  section. 

On  comparing  the  oxidation  properties  of  the  substrate  with  laser  clad 
samples.  It  was  found  that  all  laser  clad  samples  show  better  oxidation 
properties  and  especially  under  focused  laser  clad  samples  at  oxidation 
temperature  1050°C.  In  order  to  understand  the  high  temperature  oxidation 
behavior  of  these  alloys,  microstructural  investigation  of  the  laser  clad 
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samples  were  carried  out.  The  optical  micrographs  of  the  laser  clad  sample 
are  presented  in  Fig.  4.  It  shows  uniform  distribution  and  formation  of 
second  phase  in  the  laser  clad  matrix.  As  there  was  not  any  sharp  interface 
boundary  between  the  laser  clad  region  and  substrate  matrix,  it  suggests  a 
good  mixing  and  strong  metallurgical  bonding  between  the  laser  clad  region  and 
substrate.  Figures  5a  and  b  are  the  optical  and  SEM  micrographs  of  the  laser 
clad  sample  (at  5  kW,  6.35  mm/sec  with  over  focused  laser  beam)  reveals  the 
uniform  distribution  of  undissolved  Hf  and  Hf  rich  precipitates  throughout  the 
matrix.  The  segregation  of  undissolved  Hf  particles  were  not  observed  at  the 
grain  boundaries.  By  changing  the  laser  processing  conditions  (i.e.  laser 
power  6  kW  and  traverse  speed  6.5  mm/sec  with  under  focused  laser  beam),  the 
size  and  volume  of  second  phase  precipitate  were  also  changed.  Fig.  5c  and 
d.  Relatively  large  volume  fraction  with  fine  distribution  of  second-phase 
precipitate  as  well  as  undissolved  Hf  particles  were  observed.  This  is  due  to 
mainly  laser  optics  interaction  with  the  substrate  during  melting  which  would 
be  discussed  in  the  next  section.  To  see  the  detailed  internal  structure  of 
the  laser  clad  region,  STEM  investigation  was  carried  out.  The  general  survey 
of  the  thin  foils  show  fine  and  uniform  distribution  of  undissolved  Hf 
particle  as  marked  by  arrows,  and  also  Hf  rich  precipitates  in  the  nickel 
matrix  (Fig.  6a  and  b).  In  addition,  very  fine  gamma  prime  (y1)  precipitates 
were  also  observed  In  the  matrix  (Fig.  6c).  The  average  size  of  y'  precipi¬ 
tate  is  about  -70  A°.  The  y'  precipitates  are  coherent  with  the  matrix  and 
can  only  be  seen  In  the  dark  field  micrograph  using  one  of  the  y'  precipitate 
super  lattice  reflections  (Fig.  6c).  The  average  composition  of  the  matrix  is 
displayed  in  Fig.  7.  Figure  8  displays  an  electron  micrograph  of  the  laser 
clad  region  showing  the  precipitates  containing  high  density  of  stacking 
faults,  (Fig.  8a  and  b)  and  corresponding  diffraction  pattern  from  the  precip- 
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itate  is  presented  in  Fig.  8c.  The  average  composition  of  precipitate  is 
presented  in  Fig.  9,  which  reveals  that  the  precipitate  is  rich  in  Hf.  As  the 
equilibrium  phase  diagram  of  the  Ni-Cr-Al-Hf  quarternary  system  is  not 
available  and  often  metastable  phases  are  formed  during  rapid  solidification 
process  [11],  so  it  is  difficult  to  comment  about  the  exact  crystal  structure 
of  this  Hf  rich  phase.  As  the  precipitates  contains  mainly  Hf  and  Ni  elements 
with  very  small  amount  of  Cr  and  A1 ,  thus  pseudo  binary  Ni -Hf  diagram  can  be 
considered.  Based  on  the  diffraction  analysis  (Fig.  8c)  and  microchemistry  of 
the  precipitate  (Fig.  9),  it  appeared  that  the  precipitate  is  having  hexagonal 
crystal  structure  [11]. 


Discussion 

It  has  been  observed  that  inherent  rapid  melting  and  solidification 
during  laser  cladding  process  produces  very  fine  grained  microstructure  with 
increased  the  solid  solubility  of  Hf  in  matrix.  The  clad  alloys  were  found  to 
have  better  oxidation  resistance  properties  as  compared  with  the  substrate  and 
this  is  due  to  the  combined  effect  of  increased  solid  solubility  of  Hf,  Hf 
rich  precipitate  and  undissolved  Hf  in  the  Ni  matrix  (Fig.  3  and  Table  I).  At 
this  stage,  it  is  difficult  to  comment  about  which  particular  mechanism  or 
groups  of  mechanisms  are  responsible  for  improving  high  temperature  oxidation 
resistance  properties,  as  various  mechanisms  were  proposed  by  various  authors 
for  the  MCrAIRE  alloys  [3-8].  The  most  beneficial  effect  of  laser  cladding 
process  was  found  to  produce  fine  microstructure,  uniform  distribution  of  Hf 
rich  precipitate  and  undissolved  Hf  particles.  In  addition,  the  undissolved 
Hf  was  in  the  form  of  hafnium  and  not  the  oxide  of  hafnium  i.e.,  hafnia  [11]. 
During  oxidation,  the  undissolved  Hf  will  form  Its  oxide  which  will  probably 
act  as  a  pegs,  which  will  Increase  the  bonding  between  the  oxide  scale  and 
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substrate  [161.  Similarly,  Hf  rich  intermetall ic  precipitates  would  decompose 
at  high  temperature  and  would  form  its  oxide  during  oxidation  process.  Same 
behavior  would  be  expected  with  Hf  in  the  matrix.  Therefore,  there  would  be 
always  a  competitive  reaction  for  the  formation  of  oxide  of  hafnium  among  the 
hafnium  present  at  different  states.  It  is  believed  that  these  reactive- 
element  oxide  colonies  provide  excellent  sinks  for  excess  vacancies  in  the 
metal  lattice  that  may  have  otherwise  condensed  at  the  scale/metal  interface 
[17).  Condensation  of  vacancies  at  the  scale/metal  interface  would  promote 
spallation  of  the  scale.  Addition  of  Hf  to  the  Ni  based  super  alloys,  not 
only  improves  the  oxidation  resistance  properties,  but  it  also  improves  the 
mechanical  properties  (especially  creep  properties)  and  controls  the  growth  of 
gamma  prime  precipitates  during  high  temperature  cyclic  treatments  [181. 

Theoretical  calculations  have  shown  that,  in  nickel  based  alloys  con¬ 
taining  aluminum  with  RE  additions  (Y  or  Hf),  the  RE  atoms  have  a  strong 
binding  energy  to  the  aluminum  atoms  [191,  because  of  charge  transfer  from 
their  high-lying  valence  orbitals  to  lower-lying  Ni  s-d  bands.  Thus,  the 
outward  diffusion  of  cations  would  apparently  be  suppressed  because  Cr  atoms 
are  tied  up  by  RE  (Hf)  atoms.  Therefore,  the  mechanism  of  scale  growth  may  be 
controlled  by  the  inward  transport  of  oxygen.  Further,  as  mentioned  above, 
that  the  reactive  elements  provides  an  excellent  sites  to  condense  excess 
vacancies,  it  also  acts  as  nucleation  site  for  the  continuous  formation  of 
O2O3  layer  thus  prevent  the  depletion  of  Cr  content  in  the  alloy  adjacent  to 
the  interface  [191. 

At  this  stage,  it  is  difficult  to  comment  about  the  extent  of  contri¬ 
bution  of  the  size  and  volume  fraction  of  y'  precipitate  present  in  matrix 
towards  the  oxidation  resistance  properties.  It  has  been  observed  that  the 
volume  fraction  and  size  of  y'  precipitate  Increases  with  increasing  the  laser 
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power  or  decreasing  the  traverse  speed  120].  In  the  y'  (NigAl)  precipitate, 
A1  would  be  replaced  by  Hf  where  as  Ni  would  be  replaced  by  Cr,  Co,  Fe,  etc. 
[11].  In  order  to  have  good  binding  energy  between  Hf  with  A1 ,  Ni,  Cr,  it 
would  be  an  ideal  to  have  optimum  size  and  volume  fraction  of  y'  precipitate 
present  in  the  matrix;  which  would  probably  suppress  the  outward  diffusion  of 
cations  and  increases  the  oxidation  resistance  properties. 

The  oxidation  properties  of  under  focused  laser  clad  alloys  were  found  to 
be  better  than  over  focused  laser  clad  alloys.  This  is  due  to  concentration 
difference  of  alloying  elements  present  in  a  given  volume  of  melt  pool  region 
which  would  be  discussed  in  detail  as  follows.  The  laser  beam  diameter  (3  m) 
is  the  same  as  over  focused  and  under  focused  position  (Fig  10a),  but  their 
interaction  with  the  substrate  is  different.  Figure  10a  represents  a  simple 
laser  beam  optics.  When  the  substrate  is  positioned  at  over  focused  laser 
beam  location  with  3  mm  beam  diameter  (Fig.  10b),  the  laser  beam  is  in  the 
converging  mode.  Once  the  laser  beam  is  incident  on  the  substrate,  it  would 
like  to  penetrate  inside  the  substrate  up  to  the  focal  point.  The  laser  beam 
would  melt  the  substrate  up  to  focal  point  i.e.  depth  dj  of  substrate  and  mix 
with  melted  powder  to  form  a  clad  layer.  If  the  substrate  is  located  at  the 
under  focused  laser  beam  with  3  mm  diameter  (Fig.  10c),  the  interaction  of 
laser  beam  with  substrate  would  be  always  with  diverging  beam  with  lower  power 
density.  The  depth  of  penetration  of  laser  beam  inside  the  substrate  would  be 
low  (d£)  and  would  mix  with  same  volume  of  melted  powder  to  form  a  clad 
layer.  Since  the  melt  depth  (d^)  for  over  focused  laser  beam  is  more  than  d^ 
for  under  focused  condition,  melted  volume  for  the  over  focused  laser  beam 
will  be  higher.  This  suggests  that  laser  clad  alloy  produced  by  over  focused 
laser  beam  would  be  more  diluted  with  alloying  elements  such  as  Cr,  A1  and 
Hf.  This  will  also  directly  control  the  total  volume  fraction  of  Hf  rich 
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precipitates  formed,  solid  solubility  of  Hf  in  the  matrix  and  undissolved  Hf 
particles  present  in  the  total  volume  of  solidified  melt  pool  region.  Figure 
5  supports  the  above  statement  that  total  volume  fraction  of  Hf  rich 
precipitates  and  undissolved  Hf  particles  are  more  in  case  of  under  focused 
laser  clad  alloys  as  compared  with  over  focused  laser  clad  alloys  because  they 
were  rich  in  alloying  contents  present  in  the  less  solidified  melt  pool 
region.  Therefore,  better  oxidation  properties  are  expected  for  laser  clad 
alloys  produced  by  under-focused  laser  beam. 


Summary 

By  laser  cladding  process,  fine  microstructure  with  uniform  distribution 
of  undissolved  hafnium,  Hf  rich  intermetallic  precipitates  and  matrix  with 
extended  solid  solubility  of  Hf  in  the  laser  clad  region  was  achieved.  The 
overall  oxidation  resistance  properties  of  the  laser  clad  alloy  produced  by 
under-focused  laser  beam  was  found  to  be  better  than  the  alloy  produced  by 
over-focused  laser  beam  as  well  as  Inconel  718  substrate.  The  present  results 
suggest  that  the  reactive  elements  would  change  the  mechanism  of  the  scale 
formation  as  a  competitive  reaction  is  taking  place  for  the  oxidation.  These 
reactive  oxidized  particles  will  act  as  sinks  for  excess  vacancies  thus 
inhibiting  their  condensation  at  the  scale/matrix  interface  and  possibly 
enhancing  scale  adhesion. 
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Table  I  Oxidation  Data  for  Inconel  718 
(60  minutes  at  900*0) 


No.  of  OBS  Initial  Wt.,  grams  Final  Wt..  grams  Wt.  Change, 


2.2856 

0.6618 


2.2843 

0.6594 


-1.3  x  kH 
-2.4  x  10"J 


Figure  1  -  Optics  used  for  laser  cladding. 
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Figure  4  -  Optical  micrograph  shows  uniform  microstructure  in 
the  laser  clad  region,  good  interaction  between 
the  laser  clad  substrate  interface. 
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Figure  5  -  Micrographs  of  laser  clad  regions  for  over-focused  (a)  and  (b)  and  under-focused 
(b)  and  (c)  laser  beam  (a)  Optical:  showing  uniform  distribution  of  second  phase 

(b)  SEM:  showing  distribution  of  Hf  rich  precipitates  and  undissolved  Hf  particles 

(c)  Optical  showing  relatively  large  volume  fraction  of  second  phase  (d)  SEM:  showing 
fine  and  large  volume  fraction  of  second  phase  as  well  as  undissolved  Hf  particles. 
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Figure  8  -  Transmission  electron  micrograph  of  the  laser  clad  sample 
showing  the  meta  stable  phases  (a)  bright  field 
(b)  dark  field  (c)  corresponding  diffraction  pattern. 
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ONE-DIMENSIONAL  DIFFUSION  MODEL  FOR  EXTENDED  SOLID 
SOLUTION  IN  LASER  CLADDING 
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ABSTRACT 

A  mathematical  model  is  presented  for  determining  the  composition  of  ex¬ 
tended  solid  solution  formed  due  to  rapid  cooling  in  laser  cladding.  This 
model  considers  diffusion  mechanism  for  mass  transport  In  a  one-dimensional 
semi-infinite  molten  pool  of  the  cladding  material  from  which  heat  is  removed 
by  conduction  through  a  one-dimensional  semi-infinite  solid  substrate.  The 
rate  of  solidification  has  been  obtained  by  modeling  the  cooling  process  as  a 
composite  medium  heat  transfer  problem.  The  discontinuity  of  the  concentra¬ 
tion  field  has  been  simulated  using  a  nonequilibrium  partition  coefficient  and 
then  a  non-similar  exact  solution  for  the  mass  transport  equation  has  been  ob¬ 
tained  using  a  set  of  similarity  variables  which  has  been  derived  using  Lie 
group  theory. 


I.  INTRODUCTION 


Inherent  rapid  heating  and  cooling  rate  in  laser  surface  modification 
provides  an  opportunity  to  produce  novel  materials  without  being  restricted  by 
equilibrium  phase  diagram.  In  laser  cladding  and  alloying  processes  surface 
chemistry  of  materials  can  be  significantly  altered  by  simultaneous  addition 
of  alloying  elements  and  their  subsequent  rapid  solidification.  High  cooling 
rate  (=  10^°K/sec)  often  leads  to  extension  of  solubility  of  the  solute  atoms 
and  thus  produce  novel  metastable  materials.  This  process  is  dependent  on 
energy,  momentum  and  mass  transport.  Energy  transport  determines  the  rate  of 
heating  and  cooling  whereas  momentum  and  mass  transport  determine  the  extent 
of  mixing  and  final  composition. 

This  paper  is  concerned  with  mathematical  modeling  of  solute  redistribu¬ 
tion  during  laser  cladding.  The  cladding  material  which  could  be  a  single  el¬ 
ement  or  a  mixture  of  several  elements  is  delivered  to  the  surface  of  a  solid 
substrate  moving  at  a  constant  velocity,  and  is  melted  simultaneously  using  a 
laser  beam.  The  molten  pool  of  cladding  material  cools  down  mostly  by  con¬ 
ducting  away  heat  through  the  solid  substrate  and  by  losing  some  of  its  energy 
to  the  ambient  air  by  radiation  and  convection.  A  very  good  heat  transfer 
model  for  CW  laser  material  processing  was  developed  by  Mazumder  and  Steen*. 
Convection  arising  due  to  surface  tension  in  a  laser  melted  pool  was  con¬ 
sidered  by  Chan,  et  al.2.  They  studied  its  effect  on  surface  velocity,  sur¬ 
face  temperature,  pool  shape  and  the  cooling  rate.  The  distribution  of  matte'* 
by  diffusion  and  convection,  after  It  Is  delivered  to  such  a  pool  was  examined 
by  Chande  and  Mazumder^. 

There  are  not  very  many  theoretical  studies  published  In  the  literature 
to  model  laser  surface  cladding  and  alloying  at  a  nonequil ibrium  cooling 
rate.  For  equilibrium  partitioning  of  solute  at  the  solid-liquid  Interface  of 


a  binary  alloy  Tiller4  et  al.  provided  a  quantitative  analysis  for  transient 
and  steady  state  conditions.  In  this  work  they  avoided  the  heat  transfer 
analysis  by  assuming  that  the  process  of  solidification  is  taking  place  at  a 
constant  speed.  Tien  and  Geiger^*®  performed  a  heat  transfer  analysis  of  the 
solidification  of  a  semi-infinite  binary  eutectic  system.  Muelhbauer7  et  al. 
obtained  approximate  solutions  for  the  temperature  distribution  and  the  rate 
of  phase  change  for  one-dimensional  solidification  of  a  finite  slab  of  a  bi¬ 
nary  alloy.  They  considered  a  non-eutectic  system  where  solidification  takes 

Q 

place  over  a  range  of  temperatures.  Tsubaki  and  Boley  carried  out  both  heat 
and  mass  transfer  analysis  for  binary  mixtures  with  concentration  dependent 
freezing  point  and  equilibrium  partitioning  of  solute  at  the  solid-liquid  in¬ 
terface.  In  all  these  studies  the  objectives  were  either  to  analyze  only  the 
process  of  solidification  under  equilibrium  as  well  as  nonequilibrium  cooling 
rate  or  to  study  the  heat  and  mass  transfer  under  equilibrium  conditions  using 
space-time  variables.  There  has  also  been  a  considerable  effort  to  predict 
the  formation  of  extended  solid  solution  under  rapid  cooling  using  thermo¬ 
dynamic  variables  such  as  free  energy  and  chemical  potential.  The  thermo¬ 
dynamic  approach  to  solidification  has  been  examined  very  well  by  Baker  and 
Cahn9.  Boettinger  and  Perepezko10  have  discussed  the  process  of  rapid  solidi¬ 
fication  from  the  point  of  view  of  thermodynamics.  Boettinger  et  al.11  have 
also  used  the  response  function  approach  of  Baker  and  Cahn*  and  stability 
analysis  for  ml crosegregat Ion-free  solidification.  More  about  rapidly 
solidified  materials  can  be  found  In  Ref.*2’*2.  The  present  study  is  con¬ 
cerned  with  the  determination  of  composition  of  the  extended  solid  solution 
formed  due  to  rapid  cooling  in  laser  cladding  using  space-time  variables. 


II.  MATHEMATICAL  MODEL 
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Process  Physics 

The  physical  process  of  laser  cladding  involves  pouring  of  cladding 
material  powder  onto  the  substrate  to  be  clad  by  using  a  powder  delivery  sys¬ 
tem  and  simultaneously  melting  the  powder  using  a  laser  beam.  The  pool  of  the 
cladding  melt  which  is  formed  just  below  the  laser  beam  solidifies  by  losing 
heat  to  the  surrounding  air,  adjacent  cladding  and  the  solid  substrate  as  it 
moves  away  from  the  laser  beam  due  to  the  motion  of  the  substrate.  This  pro¬ 
duces  a  layer  of  cladding  on  the  solid  substrate  and  a  liquid  cladding  pool 
below  the  laser  beam  (see  Fig.  1(a)).  Figure  1(a)  shows  an  experimental  set-up 
for  laser  surface  cladding.  The  area  ABCA  shows  the  liquid  pool  of  cladding 
melt  which  is  in  contact  with  the  solid  substrate  along  AB.  It  is  also  in 
contact  along  BC  with  the  cladding  formed  on  the  substrate  and  the  portion  CA 
of  the  pool  Is  exposed  to  an  inert  atmosphere  to  avoid  oxidation  of  the 
cladding  material.  The  shape  of  the  strip  of  cladding  will  depend  on  the 
properties  and  the  temperature  of  the  cladding  melt  and  the  substrate, 
relative  motion  of  the  cladding  with  respect  to  the  substrate  and  the 
characteristics  of  the  cladding  powder  delivery  system.  In  the  present  model 
the  strip  of  cladding  BCDEB  (see  Fig.  1(a))  has  been  considered  to  have  semi- 
cylindrical  shape  (see  Fig.  1(b))  and  based  on  this  the  cladding  radius  (rc) 
is  determined  (see  Eq.  (18)).  This  paper  is  concerned  with  heat  and  mass 
transfer  In  the  cladding  melt  ABCA  to  model  extended  solid  solution  during 
laser  cladding.  This  liquid  pool  solidifies  by  conducting  away  heat  to  the 
substrate  and  to  the  solid  cladding  across  AB  and  BC,  respectively.  Also,  it 
loses  some  energy  to  the  ambient  Inert  gas  across  the  free  surface  CA.  But 
the  loss  of  energy  to  the  Inert  gas  across  the  surface  CA  will  be  smaller  than 
the  heat  loss  due  to  conduction  in  the  substrate  and  the  cladding  BCDEB. 
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However,  the  cladding  BCDEB  is  formed  from  the  liquid  cladding  melt  just 
solidified  and  hence  its  temperature  distribution  is  almost  uniform  and  much 
higher  than  the  temperature  of  the  substrate  which  is  immediately  below  the 
pool  ABCA.  So  the  heat  loss  across  BC  is  much  less  compared  to  that  across 
AB.  Therefore,  we  can  consider  that  the  pool  of  cladding  melt  (ABCA)  loses 
energy  only  in  one  direction  through  the  surface  at  AB  and  to  simplify  the 
heat  transfer  calculations  we  assume  that  the  pool  extends  up  to  infinity 
along  AB  as  well  as  BC.  Due  to  this  assumption  the  cladding  melt  solidifies 
along  AB  and  this  freezing  front  moves  upward  in  the  direction  of  BC.  The 
solid  substrate  has  also  been  taken  to  be  semi-infinite.  The  geometric 
configuration  of  this  pool  has  been  shown  in  Fig.  1(c)  by  rotating  the  pool  90 
degrees  in  the  clockwise  direction.  In  this  figure  the  freezing  front  has 
been  shown  to  be  planar.  This  is  true  when  a  pure  metal  solidifies.  For  a 
cladding  melt  of  an  alloy  system,  the  freezing  front  develops  curvature  due  to 
surface  tension  and  also  there  could  be  dendrite  and  cellular  growth  at  the 
solid-liquid  interface.  Moreover,  rapid  quenching  may  lead  to  growth  due  to 
spinodal  decomposition  and  nucleation  in  the  bulk  of  the  liquid  phase.  The 
stability  criterion  of  Mullins  and  Sekerka^®  of  a  planar  interface  during 
solidification  of  a  dilute  binary  alloy  shows  that  the  planar  interface  will 
be  unstable  for  both  the  nickel-hafnium  and  the  nickel-aluminum  systems 
considered  in  this  study.  The  growth  during  solidification  for  these  two 
cladding  alloys  has  been  shown  in  Figs.  10  and  11.  Figure  10  shows  the 
solidification  of  nickel-hafnium  cladding  powder  on  nickel  substrate  and  the 
Fig.  11  shows  the  growth  of  nickel-aluminum  cladding  powder  on  nickel  sub¬ 
strate.  It  can  be  seen  from  these  two  figures  that  there  are  dendritic 
growths  for  both  cases,  although  the  dendrites  are  less  prominant  for  the 
nickel-aluminum  alloy  compared  to  those  of  nickel-hafnium  system.  Thus  the 
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solid-liquid  interface  for  nickel-aluminum  system  can  be  considered  planar. 
Also  for  those  systems  where  the  size  of  dendrites  is  extremely  small,  the 
interface  can  be  taken  to  be  planar.  Nevertheless,  the  deidrites  will  affect 
the  diffusion  of  solute  atoms  at  the  freezing  front.  To  simplify  the  mass 
transfer  analysis,  the  freezing  front  has  been  assumed  to  be  planar  because 
the  dendrite  size  is  negligible  for  systems  such  as  nickel-aluminum.  More¬ 
over,  the  trend  predicted  by  the  model  is  more  important  than  the  exact 
numerical  values  due  to  the  paucity  of  high  temperature  materials  data  and 
simplifying  assumptions.  Besides  this  the  present  model  has  been  developed 
based  on  a  few  more  assumptions  listed  below. 

(1)  The  thermal  conductivity  and  the  thermal  diffusivity  for  a  mixture  is 
the  sum  of  the  volume-averaged  value  of  the  respective  transport 
properties  of  each  element  of  the  mixture, 

(2)  The  mass  diffusivity  of  each  element  in  the  liquid  phase  is  the  aver¬ 
age  value  of  self-diffusivity  over  the  room  temperature  and  the  in¬ 

itial  temperature  with  modified  activation  energy  for  the  mixture, 

(3)  There  is  no  diffusion  of  mass  in  the  solid  phase, 

(4)  The  solute  segregated  at  the  solid-liquid  interface  moves  in  the 

liquid  phase  by  diffusion  only, 

(5)  The  concentration  of  solute  in  the  liquid  alloy  is  equal  to  the 

nominal  composition  of  the  cladding  powder  mixture. 

(6)  Change  in  solute  concentration  in  the  liquid  phase  does  not  alter  its 
freezing  point  but  it  affects  the  heat  flux  across  the  solid-liquid 
interface,  and 

(7)  Only  50  percent  of  the  laser  energy  has  been  assumed  to  be  absorbed 
by  the  cladding  material.  Studies^  show  that  the  amount  of  laser 
energy  absorbed  by  different  materials  is  37-60  percent. 
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Under  these  assumptions  it  is  required  to  solve  the  one-dimensional  heat 
conduction  equation  to  obtain  the  speed  of  the  freezing  front  which  is  then 
utilized  to  solve  for  the  distribution  of  solute  in  the  liquid  phase.  In 
laser  processing  of  material  usually  the  laser  beam  or  the  workpiece  moves  at 
a  certain  speed  and  when  the  laser  melted  pool  is  fully  developed  the  speed  of 
the  freezing  front  and  that  of  the  laser  beam  or  the  workpiece  are  identi¬ 
cal.  This  eliminates  the  need  for  solving  the  energy  transport  equation  to 
determine  the  speed  of  solidification  for  fully  developed  pool.  In  the  case 
of  laser  surface  alloying  the  solidification  progresses  with  segregation  of 
solute  atoms  from  the  solid  to  the  liquid  phase  and  this  alters  the  concentra¬ 
tion  of  solute  at  the  sol  id- liquid  interface.  This  will  affect  the  speed  of 
the  freezing  front  as  well  as  the  freezing  point  of  the  liquid.  To  take  the 
former  effect  into  account  the  energy  transport  equation  has  been  solved  while 
the  latter  one  has  been  ignored  in  this  model. 


Mathematical  Formulations 

The  governing  equations  for  energy  transport  are: 


(i)  Region  1  (: 

S 

.  i 

V. 

3X2  °1 

<  x  <  0 


(ii)Region  2  (solidified  cladding) 


3  ^2  1  3^2 

“T  =  a  t  >  0,  0  <  x  <  S(t) 

3xZ  a2  3t 


(iii)Region  3  (liquid  cladding) 


I  V*  •j'j*  M  •  m 


j.  v.v  .v ~ /• 


E! 


91 


i 

s 


3 


— r  =  TrTT'  t  >  0,  S(t)  <  X  <  -  (3) 

Here  Tj,  T2,  and  T3  are  the  temperatures  and  c^,  c^,  and  are  thermal 
diffusivities  of  regions  1,  2,  and  3,  respectively.  S(t)  is  the  location  of 

the  solid-liquid  interface  at  time  t.  If  k ^ ,  k2,  and  k3  are  taken  to  be 

thermal  conductivies  of  the  regions  I,  2,  and  3,  respectively  and  if  the 
ambient  temperature,  the  freezing  point  of  the  cladding  melt  and  the  initial 
cladding  pool  mean  temperature  are  denoted  by  T  ,  Tm,  and  Tf,  respective^ 

then  the  auxiliary  conditions  for  the  above  heat  transfer  problem  can  be 

written  as 


Tl<--  l>  ■  Tc 

(4a) 

I 

V 

N* 

T^x.O)  =  Tq 

(4b) 

s 

aT.  k_  sT_ 

— L  =  ~  — -  at  x  =  0 

3 X  k^  3X 

(4c) 

r, 

>. 

& 

T  L (0 , t)  =  T2(0,t) 

(4d) 

i7 

v, 

T2(x,t)  =  T3(x,t)  =  lm  at  x  =  S(t) 

( 4e 

i 

3T«  3T-  jc 

k2  73T-  k3TF  =  pLIat  *  -  s^) 

:  J* 

ft 


where  p  is  the  density  of  the  cladding  meU  a-vi 
solidification.  Note  that  the  change  in  jpnS  » ,  1... 
been  neglected  here. 
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T,(-.t)  =  T- 


T3(x,0)  =  Tf 


S(0)  =  0 


The  amount  of  solute  rejected  by  the  solidifying  cladding  and  its  transport  in 
the  liquid  phase  are  governed  by  the  following  diffusion  equation. 


aC  _  ^3  aC  _  p  a2C 
3t  /t  3X  ax2 


where  c  is  the  concentration  of  solute  in  the  liquid  phase,  0  is  the  mass 
diffusion  coefficient  and  y  Is  a  root  of  the  transcendental  Eq.  (12).  If  CQ 
is  the  initial  concentration  of  solute  in  the  liquid  phase  and  kg  denotes  the 
equilibrium  partition  coefficient  at  the  solid-liquid  interface  then  the 
auxiliary  conditions  for  the  above  mass  transfer  problem  are 


C(x,o)  =  C 


C(-.t)  =  C 


y/aT  a  C  */ct3 

D(1  +  -_i-  B)  ik  =  —A.  (k  -  1)C  at  x  -  0 
✓t  3  /t  e 


in  a  coordinate  frame  of  reference  moving  with  the  solid-liquid  interface.  In 
the  expression  (6c),  e  *  x/0*,  where  x  is  the  interatomic  spacing  and  D*  is 


the  inter-diffusivity.  This  expression  is  obtained  by  incorporating  into  the 
mass  balance  equation  at  the  solid-liquid  interface  an  expression  for  non¬ 
equilibrium  partition  coefficient  derived  by  Aziz*^. 

Method  of  Solution 

The  multi-region  heat  transfer  problem  can  be  solved  by  the  method  used 
by  Carslaw  and  Jaeger*4.  The  temperature  distributions  in  various  regions  are 
given  by 


V"1> 

^2(^2* 

T3<"3> 


“iMt„  -  T„> 

To  *  a?kjl  :  irfTVTT  11  *  Erf<"l>' 


«lk2<T*  -  V 


°2kl 


To  +  agkjl  +  Erf (y)T  11  +  Erf  (n2)] 


T_  -  T 


■  Tf  *  ErfcJVJ  Erfc<"3> 


(7) 

(8) 

(9) 


where 


n..  =  x/2/a.t,  i  =  1,2,3 

and  the  location  and  the  speed  of  the  freezing  front  are  given  by 


set)  =  2,  (.3t)1/2 

(10) 

S(t)  =  y/a^/Zt 

(11) 

Here  k£  and  are  taken  to  be  equal  to  k^  and  respectively  and  y  is  ob¬ 
tained  from  the  following  transcendental  equation: 
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After  knowing  the  rate  of  formation  of  the  solid  phase  the  distribution  of 
solute  in  the  liquid  phase  can  be  obtained  by  solving  the  mass  diffusion  prob¬ 
lem  described  by  Eq.  (5)  and  the  associated  conditions  (6a-c).  An  exact  solu¬ 
tion  for  the  above  problem  can  be  obtained  by  using  similarity  variables  which 
are  determined  using  Lie  group  theory.  Details  of  this  theory  and  its 
application  for  the  solution  of  differential  equations  can  be  found Us¬ 
ing  this  group  theoretic  approach  the  similarity  variables  can  be  found  to  be 
(See  APPENDIX  A). 
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di|>  y/a-5  (ke-l)  y/o3  (k  -1)  n-1  i 

W  -  D6  *n  ■  f1  '  *»>  — D8  *  (-1'  J  *J  <16b> 


where  6  is  the  Kronecker  delta  and  A's  of  expression  (14)  are  given  by 
no  n 


An  =  (er/a^)  ,  n  =  0,  1,  2,  . .  ® 


For  this  problem  four  terms  of  the  series  (14)  were  computed  by  determin¬ 


ing  ii»0,  4)^,  and  from  Eq.  (15).  These  solutions  are  (see  APPENDIX  B) 


ib  (r)  =  C  [1 - — - G  Erfc  _(£ 2  r) - - - 

0  0  (/8/ir)  Exp(-4/  c*3/D)  +  G  Erfc(2Y/^7D] 


Mr) 


G  *0(r0)  -<r2  -  r/)/2 


ro  +  G 


,  .  -  *l<ro>'  -<r2  -  ro2>/2 

*2(r)  * - — - - H -  •  r  e 

2  roG  +  <ro  ' 


GU  (r  )  -  Mr  )  +  4>2(r  ) ]  ~  -(r2  -  r  2)/2 

Mr)  -  -  — f-* - — - -  1)  e 

3  (ro  -  DG  +  r0<3  -  ro  > 


where 


G  =  y/a3  (kg  -  1)  /2/D,  r  =  (y  +  2Y/a3)//2D  and  rQ  =  Y/2a3 


/2a, /D 
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Using  these  expressions  for  4>q,  ^  and  ^  the  concentration  of  solute  In 
the  liquid  phase  was  computed  using  the  Eq.  (14)  and  then  the  corresponding 
concentration  In  the  solid  phase  (Cs)  was  obtained  using  the  following 
expression  derived  by  Aziz*®  for  nonequilibrium  partition  coefficient. 


S(t)s  +  kc 


S(t) 6  +  1 


In  the  present  model  the  inter-diffusivity  D*  is  taken  to  be  equal  to  the 


mass  diffusivity  and  the  interatomic  spacing  x  is  taken  to  be  4A  which  is  the 


average  diameter  of  an  atom.  Other  physical  properties  for  the  alloy  were 


computed  at  its  nominal  composition  using  the  following  relationships: 


p  =  1/  r  W./p. 


where  n  *  total  number  of  elements  present  in  the  liquid  pool  of  cladding 


L  =  r  W.  L. 
j=l  J  J 


n  W 

k  (T)  -  O  £  -1  k..(T) 
1  j=l  pj 


k1  =  J  k1(T)dT/(Tf  -  Tq) 
o 


^(T)  =  r  W.  Z  (T) 


'r«v 1  /  "1  ''n  ' 

P  >1  J  Pj 


Cp  ■  {  ^pCT)dT/(Tm  -  Tq) 
0 


n  W.  _ 

a.(T)  =  p  £  a.^T) 

1  J-l  °j  1J 


S 


B 


B 
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'f  __ 

a.  =  J  oi(T)dT/(Tf  -  Tq) 

0 

and  the  diffusion  coefficient  for  the  j-th  solute  in  the  liquid  phase  is 

Tf  -  L 

°J  *  {  °oj  e  RT  dT/<Tf  -  To> 

0 

where  DQj  is  the  self-diffusion  coefficient  of  the  j-th  element  and  Q  is  the 
activation  energy  for  the  mixture  given  by 

n 

Q  =  t  w.Q. 
j-1  J  J 

In  all  these  expressions  for  thermophysical  properties  of  the  cladding  powder 
mixture,  the  index  i  =  1,2,3  refers  to  the  solid  substrate,  solidified 
cladding  and  the  liquid  cladding  region,  respectively,  and  j  *  1,2,3,..., 

refers  to  the  j-th  element  of  the  cladding  powder  mixture.  Here,  Wj  and  wj 
are  the  weight  fraction  and  the  mole  fraction  of  the  j-th  element,  respec¬ 
tively.  k^j  and  denote  respectively  the  thermal  conductivity  and  the 

thermal  dlffusivlty  of  the  j-th  element  in  the  i-th  region.  C  is  the 

Pj 

specific  heat  of  the  j-th  element.  The  thickness  of  the  cladding  can  be 

determined  by  taking  an  over-all  energy  balance  which  yields  the  following 
expression  for  the  radius  (rQ)  of  the  semi-cylindricel  strip  of  cladding: 


where  f  is  the  fraction  of  laser  energy  absorbed  by  the  cladding  material,  P 
is  the  power  of  the  laser  beam  Incident  on  the  material,  and  v  is  the  speed  of 
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the  workpiece  which  can  be  related  to  the  volumetric  rate  (q)  of  powder 
delivery  by 


q 


2 

*rCv 


(19) 


assuming  that  all  of  the  delivered  powder  form  cladding  on  the  substrate. 


III.  RESULTS  AND  DISCUSSION 
Experimental  Verification: 

The  results  of  the  above  model  was  compared  with  experimental  data. 
Laser  cladding  was  performed  on  nickel  substrate  with  two  different  cladding 
materials  of  nominal  composition  74 %  Ni  and  26%  Hf  in  one  case  and  74%  Ni  and 
26%  A1  in  the  other  case.  The  solid  substrate  was  clad  with  these  cladding 
materials  which  take  the  shape  of  an  approximately  semi-cylindrical  strip  of 
metal  on  the  substrate  (see  Fig.  1(b)).  Electron  Probe  Micro  Analysis  of 
these  samples  shows  the  concentration  of  Hf  and  A1  in  the  matrix  of  Ni  in 
excess  of  what  is  predicted  by  the  equilibrium  phase  diagram.  The  present 
model  predicts  the  composition  of  the  extended  solid  solution  quite  well. 
These  results  have  been  presented  in  the  Table  I  for  Tf  =  2035°K  for  the  N i -Hf 
alloy  and  Tf  =  I907°K  for  Ni-Al  alloy. 


Table  I  Comparison  of  Theoretical  Results  with  Experimental 
Results  of  Extended  Solid  Solution 


Laser  Laser  Beam  Speed  of 
Power  Diameter  the  Work- 
(Kw)  (mm)  piece 

( Inch/mi n) 


Nominal  Composition 
of  the  cladding 
material 
(wt.  %) 

74%  Ni ,  26%  Hf 
74%  Ni,  26%  A1 


Composition  in  the  solid 
solution  (wt.  %) 

Theoretical  Experimental 
Results  Results 


3.05%  Hf 
27.3%  A1 


3.58%  Hf 
22%  A1 
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The  prediction  of  composition  of  the  extended  solid  solution  based  on  the 
present  model  is  within  30%  of  the  experimental  results.  This  discrepancy  may 
be  due  to  some  of  the  assumptions  of  this  model.  For  example,  solidification 
has  been  assumed  to  take  place  at  a  constant  temperature  which  has  been  taken 
to  be  the  freezing  point  of  the  cladding  material  at  its  nominal  compos¬ 
ition.  But  it  is  well  known  that  the  change  in  composition  in  the  liquid 
phase  will  alter  its  freezing  point.  Also  the  present  model  utilizes  an 
expression  (17)  for  nonequilibrium  partition  coefficient  which  is  applicable 
to  dilute  solution.  Apart  from  these  the  presence  of  a  two-phase  zone  between 
the  solidus  and  the  liquidus  lines  and  the  surface  tension  driven  flow  causing 
convection  in  the  liquid  pool  will  affect  the  mixing  of  solute  in  the  liquid 
phase  and  thus  alter  its  composition  in  the  solid  phase.  Paucity  of  high 
temperature  liquid  metal  data  also  contributes  to  the  numerical  error. 
However,  the  objective  of  this  study  is  to  determine  the  trend  and  thus 
understand  the  underlying  process  physics. 

Parametric  Results 

The  above  model  was  used  to  study  the  effect  of  various  important  process 
parameters  such  as  the  cooling  rate,  cladding  powder  delivery  rate  on  the  com¬ 
position  of  hafnium  in  nickel  matrix.  Results  were  obtained  for  ambient  tem¬ 
perature  Tq  =  293°K  and  some  of  the  process  parameters  were  computed  based  on 
a  laser  beam  of  diameter  3  mm.  Besides  this  the  initial  temperature  of  the 
cladding  melt  was  chosen  as  an  independent  parameter  for  one  set  of  results 
and  the  cladding  thickness  for  another  set  of  results  for  different  values  of 
laser  power  and  the  velocity  of  the  workpiece.  These  parameters  have  been 
listed  in  Tables  I  and  II. 


Table 

II  List  of  independent  parameters  (v, 
and  dependent  parameters  (rc,  q) 

P,  Tf) 

Independent 

Parameters 

Workpiece 

Velocity 

(Inch/min) 

30 

Laser 

Power 

(kW) 

5 

6 

7 

Initial 

Temperature 

(°K) 

2000  2200 

2400 

2000 

2200 

2400 

2000 

2200 

2400 

Calculated 

Values 

Cladding 

Thickness 

(mm) 

3.56  3.41 

3.28 

3.90 

3.74 

3.59 

4.22 

4.04 

3.88 

Powder 

Delivery  Rate 
(mnr/sec) 

252.8  232.0 

214.6 

303.4 

279.0 

257.1 

355.3 

325.6 

300.3 

Independent 

Parameters 

Workpiece 

Velocity 

(Inch/min) 

40 

Laser 

Power 

(kW) 

5 

6 

7 

Initial 

Temperature 

(°K) 

2000  2200 

2400 

2000 

2200 

2400 

2000 

2200 

2400 

Calculated 

Values 

Cladding 

Thickness 

(mm) 

3.20  2.96 

2.84 

3.38 

3.24 

3.11 

3.65 

3.50 

3.36 

Powder 

Oelivery  Rate 
(mnr/sec) 

255.6  233.1 

214.5 

303.9 

279.2 

257.3 

354.4 

325.8 

300.3 
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Table  III  List  of  independent  parameters  (v,  p,  rc)  and  dependent  parameters  (Tf ,q) 


Independent  Workpiece 
Parameters  Velocity 
( Inch/mi n) 

Laser 

Power 

(kW) 


Cladding 

Thickness  3.45  3.35  3.25  3.95  3.80  3.65  4.15  4.00  3.85 

(mm) 


Calculated  Initial 

Values  Temperature  2149  2292  2449 
(°K) 


Powder 

Delivery  Rate  237.4  223.9  210.7 
(mnr/sec) 


Independent  Workpiece 
Parameters  Velocity 
( Inch/min) 

Laser 

Power 

(kW) 


1948  2123  2320  2072  2247  2443 


311.3  288.1  265.8  343.6  319.2  295.7 


Cladding 

Thickness 

(mm) 


3.00  2.85  3.15  3.45  3.30  3.15  3.70  3.55  3.40 


Calculated  Initial 

Values  Temperature  2129  2383  1911 

(°K) 


Powder 

Delivery  Rate  263.9  239.4  216.1 
(mnr/sec) 


1911  2110  2337  1942  2129  2341 


316.6  289.7  263.9  364.1  335.2  307.5 
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Table  III  (Continued) 


Independent  Workpiece 
Parameters  Velocity 
(Inch/min) 


Independent 

Parameters 


Laser 

Power 

m 


Cladding 

Thickness 

(mm) 


2.75  2.65  2.60  3.00  2.90  2.80  3.30  3.15  3.00 


Calculated  Initial 

Values  Temperature  2016  2189  2382 
( °K) 


Powder 

Delivery  Rate  251.4  233.5  216.2 
(mnr/sec) 


2035  2194  2369  1954  2167  2142 


299.2  279.6  260.7  362.1  329.9  299.2 
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Results  obtained  based  on  the  parameters  of  the  Table  II  have  been  presented 
in  Figs.  2  through  4.  As  can  be  seen  in  Fig.  2  the  cooling  rate  at  the  solid- 
liquid  interface  drops  sharply  to  almost  equilibrium  cooling  rate  within  one 
second  after  the  cladding  melt  starts  freezing.  Also  it  can  be  observed  that 
the  magnitude  of  the  cooling  rate  of  the  liquid  phase  at  the  freezing  point  Is 
higher  for  lower  initial  temperature  of  the  liquid  cladding  material.  This  is 
to  be  expected  because  when  the  initial  temperature  Tf  is  lower  the  amount  of 
thermal  energy  of  the  melt  is  less  compared  to  when  it  is  at  a  higher  temper¬ 
ature.  So  when  cladding  materials  of  two  different  temperatures  are  brought 
in  contact  with  a  solid  substrate  of  temperature  293°K,  the  melt  of  lower  tem¬ 
perature  will  have  to  lose  lesser  amount  of  energy  than  the  one  which  is  at  a 
higher  temperature.  For  example,  let  us  take  a  melt  of  temperature  2000°K  and 
another  one  of  2200°K.  When  they  are  brought  in  contact  with  a  solid  sub¬ 
strate  of  temperature  293°K  the  first  melt  will  cool  down  from  2000°K  to  the 
freezing  point  by  conducting  away  its  heat  through  the  substrate  with  an  in¬ 
itial  driving  force  for  heat  transfer  proportional  to  2000  -  293  =  I707°K. 
The  second  melt  on  the  other  hand  will  first  require  some  time  to  cool  from 
2200° K  to  2000° K  during  which  time  the  temperature  of  the  solid  substrate  will 
rise  from  its  initial  temperature  293°K.  As  a  result  of  this  the  melt  will 
now  have  to  cool  form  2000°K  to  the  freezing  point  with  an  initial  driving 
force  smaller  than  that  for  the  first  melt. 

The  extension  of  solubility  In  solid  solution  has  also  been  studied  as  a 
function  of  cooling  rate.  It  can  be  seen  from  Fig.  3  that  the  higher  the 
cooling  rate  the  higher  Is  the  composition  of  hafnium  in  nickel  matrix.  How¬ 
ever  It  should  be  noted  that  same  cooling  rate  for  different  initial  melt  tem¬ 
perature  does  not  result  In  the  same  composition  of  hafnium  in  the  solid  solu¬ 
tion.  This  is  because  the  higher  the  initial  temperature  the  longer  it  takes 


to  first  form  the  solid  phase  at  the  interface  of  the  substrate  and  the  liquid 
cladding  material  and  during  this  time  the  solute  atoms  diffuse  away  from  the 
interface  causing  a  formation  of  alloy  lean  with  solute  atoms.  This  is  true 
if  there  is  no  adsorption  between  the  solute  atoms  and  the  solid  substrate. 
Another  important  process  parameter  considered  in  this  study  is  the  powder  de¬ 
livery  rate.  Its  effect  on  the  concentration  of  hafnium  for  different  laser 
energy  fluxes  has  been  shown  In  Fig.  4.  The  shift  in  the  concentration  curves 
to  the  right  as  the  laser  energy  flux  and  the  cladding  powder  delivery  rate 
are  increased  confirms  the  fact  that  an  increasing  powder  delivery  rate  causes 
more  cladding  material  on  the  substrate  which  requires  more  laser  energy  to 
melt  in  order  to  yield  a  pre-determined  concentration  of  solute  in  the 
alloy.  Powder  delivery  rate  is  a  very  important  process  parameter  since  it  is 
related  to  the  cladding  thickness  and  the  velocity  of  the  workpiece  as  well  as 
the  initial  temperature  of  the  cladding  material.  It  is  equal  to  ^  nr£v  and 
using  this  in  the  expression  for  rc  given  above  it  can  be  related  to  Tf.  Thus 
if  the  powder  delivery  rate  is  known  then  is  also  known  and  then  the  veloc¬ 
ity  of  the  workpiece  can  be  determined  for  a  desired  cladding  thickness.  It 
can  also  be  related  to  laser-cladding  interaction  time  by  noting 
that  j*r£v  =  ^itr£B/(B/v)  which  is  a  ratio  of  cladding  volume  based  on  the  beam 
diameter  (B)  to  the  interaction  time.  In  otherwords  the  powder  delivery  rate 
is  equal  to  the  cladding  volume  based  on  the  beam  diameter  for  unit  inter¬ 
action  time. 

Also,  the  weight-percentage  of  hafnium  in  the  nickel-hafnium  alloy  and 
the  initiated  temperature  of  the  cladding  melt  were  computed  for  a  set  of  in¬ 
dependent  parameters  (v,  P,  r^)  listed  in  Table  III.  These  results  have  been 
shown  In  Figs.  5  through  9,  for  various  values  of  laser-cladding  interaction 
time  and  laser  energy  flux.  Figure  5  provides  a  set  of  very  important  results 


that  are  useful  for  selecting  the  process  parameters  to  clad  a  material  with 
an  alloy  of  pre-determined  composition  and  thickness.  For  example,  the  choice 
of  the  cladding  thickness  and  the  composition  of  the  cladding  alloy  uniquely 
defines  a  curve  as  shown  in  Fig.  5  which  corresponds  to  a  particular  value  of 
the  laser  energy  flux  and  the  interaction  time.  The  latter  two  parameters  are 
related  to  the  laser  power,  p,  laser  beam  diameter,  B,  and  velocity  of  the 
workpiece,  v.  Thus  if  one  of  the  three  process  parameters,  P,  B,  or  v,  is  se¬ 
lected  along  with  rc  and  the  composition  of  the  cladding  alloy,  then  the  other 
two  parameters  can  be  obtained  from  Fig.  5.  After  this,  the  volumetric  rate 
of  cladding  powder  delivery  can  be  determined  from  Eq.  (19).  Figures  6,  7,  8 
and  9  show  the  effect  of  cladding  thickness  on  the  initial  temperature  of  the 
cladding  melt  and  the  effect  of  the  latter  variable  on  the  cooling  rate.  As 
can  be  seen  from  Fig.  6,  the  initial  temperature  decreases  with  the  increase 
in  rc  for  the  same  laser  energy  flux  and  the  interaction  time  and  Figs.  7,  8 
and  9  show  that  the  lower  the  initial  melt  temperature,  the  higher  the  cooling 
rate.  Thus  if  a  higher  cooling  rate  is  needed,  the  initial  melt  temperature 
has  to  be  chosen  from  the  results  of  Figs.  7,  8  and  9  and  then  Fig.  6  has  to 
be  used  to  determine  the  cladding  thickness. 


CONCLUSION 


The  present  work  examines  the  extension  of  solid  solubility  based  on  the 
transport  of  energy  and  mass.  Solute  transport  has  been  considered  to  take 
place  only  in  the  liquid  phase  while  the  energy  transport  has  been  considered 
in  both  solid  and  liquid  phases.  Effect  of  nonequilibrium  cooling  rate  on 
solute  segregation  at  the  freezing  front  has  been  taken  into  account  by  con¬ 
sidering  a  nonequilibrium  partition  coefficient.  Using  this  the  mass  transfer 
problem  has  been  solved  analytically  in  a  frame  of  reference  moving  with  the 
freezing  front  with  a  velocity  determined  from  the  heat  transfer  problem. 
These  mathematical  solutions  have  been  utilized  to  study  the  effect  of  various 
process  parameters  on  the  concentration  of  solute  in  an  alloy.  As  mentioned 
above  the  cooling  rate  is  found  to  drop  sharply  within  one  second  after  the 
solidification  starts  and  it  is  also  higher  for  lower  initial  melt  temper¬ 
ature.  The  hafnium  concentration  in  the  nickel-hafnium  alloy  is  found  to  be 
higher  for  higher  cooling  rate  and  also  more  laser  energy  is  required  to  ob¬ 
tain  an  alloy  of  a  given  composition  at  a  higher  powder  delivery  rate  than 
that  for  a  lower  rate  of  powder  delivery. 

Also  it  can  be  seen  from  Fig.  4  that  the  composition  of  hafnium  does  not 
depend  on  the  parameters  called  laser-cladding  interaction  time,  defined  by 
the  ratio  of  the  laser  beam  diameter  to  the  velocity  of  the  workpiece,  whereas 
it  does  depend  on  the  interaction  time  as  shown  in  Fig.  5.  This  is  because 
the  set  of  initial  temperature  is  the  same  (see  Table  II)  for  the  former  case 
but  it  is  different  (see  Table  III)  for  the  latter  one.  Another  interesting 
feature  that  can  be  observed  from  the  Fig.  5  is  that  same  compositon  of  haf¬ 
nium  is  obtained  in  the  alloy  for  different  cladding  thickness  by  varying  the 
laser  energy  flux  and  the  interaction  time.  This  is  because  the  initial  tem¬ 
perature  of  the  cladding  material  and  hence  the  cooling  rate  are  the  same  even 
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though  other  parameters  are  different  (see  Figs.  6, through  9).  So  it  can  be 
concluded  that  the  choice  of  the  initial  temperature  determines  the  compositon 
of  the  alloy.  Once  the  initial  temperature  of  the  cladding  melt  is  fixed, 
other  process  parameters  can  be  determined  by  selecting  any  two  of  the  three 
parameters,  r^,  v,  and  P  and  the  third  one  will  then  have  to  be  determined 
from  Eg.  (18). 
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APPENDIX  A 

The  similarity  variables  given  by  the  expression  (13)  can  be  obtained  by 
using  Lie  group  theory.  To  apply  this  theory  to  Eq.  (5)  let  us  define 

K  = 

and 

9  =  /t. 


This  transforms  Eq.  (5)  to  the  following  form 

K  -KCx-°Cxx*°  <A-1> 

where  the  subscripts  refer  to  the  partial  differentiation  of  C  with  respect  to 
the  subscript  variable. 

We  define  the  group  generator  for  (A.l)  as 

U  *  5(x,9,C)  +  n(x,e,C)  +  6(x,9,C)  1^-  (A. 2) 

Then  the  twice-extended  operator  becomes 


U" 


(A. 3) 


where 
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6e  +  5cCe 


C  C  £ 

X  0^C 


C0”e  -  C0  ”c 


6  +  5  C 

X  C  X 


CxC0"c 


C  C  n  -  C  ti  -  C  £ 
x  enc  e  'x  xx 


-Cx»<S  *  5c  c0>  -  C0.("e  +  "cCe>  *  see  +  s9cC0  *  sc  Ce.  + 


+  ^5ce  +  5cc^V  "  ^xe^o  '  ^x^oo  +  ('x9^c^9 


^xS^ce  +  eccCe^ 


c06ce(,e  +  <cCx>  '  VV*  *  VV  "  2ncC0C00 


co  <nce  *  nccC0> 


r "  =  -  C  (e  +  e  C  )  -  C  _ ( n  +  n  C  )  +  S 
*•  xxv^x  c  x'  X0'  x  c  x'  ) 


+  2Cxscx  +  Cx  6cc  +  6c  Cxx  -  V5xx  +  5xcCx>  - 


<x  Cxx  ~  Cx  ^cx  +  5cc  Cx>  "  2CxCxx^c  “  C9^nxx  +  nxcCx^ 


C9xK  +  Cxnc)  ~  ('9('XXnc  "  Wc*  +  nCC  ^X^ 


I 


u 


y"  =  -  C  ( £  +  £  C  )  -  C  (n  +  f)  C  )  +  6  +  <5  C  + 

'  XX  '*8  0'  X8'  8  C  8'  X8  XC  8 


+  sc  ^xe  +  ^x^ce  +  6cc  "  ^x  ^xe  ~  ^x^xe  +  ^xc 


2Cx  Cxe^c  "  Cx  («ce  +  5cc  Ce^  "  nx  Cee  “  C8^nxe  +  nxc  Ce^ 


—  n  (C  C  +  C  C  )  —  C  C  (n  +  r)  C  ) 
C  '  X  00  8  X8 '  X  0V  'C0  'cc  0' 


For  Eq.  (A.l)  to  admit  any  group,  it  is  required  that  it  be  invariant  relative 
to  the  twice-extended  operator,  that  is 


u''<i  ce  *  KCx  -  oCxx>  8  0 


(A.4) 


on  the  manifold  defined  by  Eq.  (A.l). 

When  the  expression  (A.4)  is  expanded  and  the  coefficients  of 
Cx *  C0,  Cx0,  C2,  Cx,  CXCT,  Cx0Cx,  C2  and  C0C2  are  equated  to  zero,  we  obtain 
the  following  set  of  overdetermined  partial  differential  equations: 


-  +  r  -  r  * *  \  +  \  <25x  -  sc>  *  9nxx *  0 


r-?',-  K<«c  -  «x>  -  K(25x  -  5c>  -  9<2\c  -  <xx>  =  0 


1  ‘c  +  !  «c  -  2V  •  0 


2K'c  *  e(scc  -  2exc>  •  0 


V 


This  system  of  group  determining  equations  were  solved  for  the  coordinate 


functions  5,  n  and  6.  The  solution  yields  the  following  four-parameter  group 


generator: 


U  =  (C3x  +  C4)  gX  +  C3e  gQ  +  (ClC  +  C2)  9C 


which  decomposes  into  the  following  four  generators: 


3  3 

X^7+  9  ae 


The  group  U4  yields  the  widely  used  similarity  variable  x/e,  i.e.  x/"t.  How¬ 
ever  this  similarity  variable  does  not  leave  the  interface  condition  (6c)  in¬ 
variant.  To  obtain  a  set  of  similarity  variables  that  will  be  admitted  by  the 
Eq.  (5)  as  well  as  by  the  initial  and  boundary  conditions  (6a-c)  we  observe 
that  a  linear  combination  of  Lie  groups  of  a  differential  equation  is  also  a 
group  for  the  differential  equation. 

So,  we  combine  U2  and  U4  to  obtain 

U5  =  x  fx  +  9  a?  '  nc  fc*  n  =  0,1,2 . " 

which  yields  x/t  and  C(/t)n/An  as  similarity  variables.  Here  Ap  is  an  arbit¬ 
rary  constant. 
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APPENDIX  B 

The  Eq.  (15)  is  a  linear  second  order  ordinary  differential  equation  with 
variable  coefficient.  Exact  analytical  solution  for  this  equation  for  differ¬ 
ent  values  of  n  can  be  obtained  as  follows: 


Letting  y  =  (y  +  2Y/a3)//2Q,  the  Eq.  (15)  can  be  written  as 


d  %  dii> 

— iT  +  r  dr  +  n  't'n  =  n  =  0,1,2,...,® 


For  n  =  0,  we  have 


d<l> 

+  r  —■  =  0 


and  its  solution  can  be  easily  found. 


For  n  =  1,  we  have 


2 

d  ip,  d\i>. 


which  can  be  written  as 


a?  <sr>  *  a? (r  *i>  * 0 


for  which  the  solution  is 


^(r)  =  [H1  +  Fx  f  e&2/2  del  e_f2/2 
ro 


(B.  1) 


(B.2) 


(B  .3) 


where  and  Fj  are  two  arbitrary  constants. 


For  n  =  2,  we  have 


2 

d  iiip  d<i i- 

^r  +  rdr  +  2^2-0 


(B.4) 


To  solve  this  equation  we  first  differentiate  (B.l)  with  respect  to  r  to  get 


^~2  +  r  Hr  +  (n+1)P  =  0 
dr6  ar 


(B.5) 


where 


P  = 


di|>r 

dr" 


Now 


the  Eqs.  (B.4)  and  (B.5)  are  identical  if  n>2  -  p  and  n  =  1.  So 


dill, 

h  =  dr" 


which  can  be  determined  from  the  solution  ^(r).  In  this  way  the  solutions  of 
(B.l)  can  be  obtained  from 
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Figure  1(a)  Schematic  diagram  of  laser  surface  cladding 

(b)  Three-dimensional  view  of  the  cladding  and  the  substrate 

(c)  Geometric  configuration  used  in  the  present  model.  (The 
model  substrate  and  the  solidification  of  the  mode  1  cladding 
have  been  shown  after  rotating  the  pool  ABCA  (see  Fig.  1(a))  by 
90  degrees  clockwise) 


Figure  2  Variation  of  cooling  rate  of  the  liquid  phase 
with  time  at  the  freezing  front 


Figure  3  Effect  of  cooling  rate  on  the  composition  of 

hafnium  in  the  solid  phase  at  the  freezing  front 


Figure  4  Parametric  dependence  of  the  composition  of  hafnium  in 
the  solid  phase  at  the  freezing  front  at  0.2  sec.  after 
solidification  begins,  on  laser  energy  flux 


Figure  5  Variation  of  hafnium  composition  at  the  freezing  front 
in  nickel-hafnium  alloy  relative  to  the  cladding 
thickness  at  0.2  sec.  after  solidification  begins 
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Figure  7  Effect  of  initial  temperature  on  the  cooling  rate 
at  0.2  sec  after  solidification  begins  for 
laser-cladding  interaction  time  0.24  sec 


Figure  8  Effect  of  initial  temperature  on  the  cooling  rate 
at  0.2  sec  after  solidification  begins  for 
laser-cladding  interaction  time  0.18  sec 


Figure  9  Effect  of  Initial  temperature  on  the  cooling  rate 
at  0.2  sec  after  solidification  begins  for 
laser-cladding  Interaction  time  0.14  sec 


Figure  10  Microstructure  of  nickel-hafnium 
cladding  on  nickel  substrate 


Figure  11  Microstructure  of  nickel-aluminum 
cladding  on  nickel  substrate 
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Figure  9  Effect  of  initial  temperature  on  the  cooling  rate 
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